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ABSTRACT 

We present a sample of 218 new quasar pairs with proper transverse separations i? prop < 1 hr 1 Mpc 
over the redshift range 0.5 < z < 3.0, discovered from an extensive follow up campaign to find com- 
panions around the Sloan Digital Sky Survey and 2dF Quasar Redshift Survey quasars. This sample 
includes 26 new binary quasars with separations R pTO p < 50 kpc (6 < 10"), more than doubling 
the number of such systems known. We define a statistical sample of binaries selected with homoge- 
neous criteria and compute its selection function, taking into account sources of incompleteness. The 
first measurement of the quasar correlation function on scales 10 h" 1 kpc < i? pr0 p < 400 hr 1 kpc is 
presented. For i? prop < 40 hr 1 kpc, we detect an order of magnitude excess clustering over the expec- 
tation from the large scale (i? prop > 3 h^ 1 Mpc) quasar correlation function, extrapolated down as a 
power law to the separations probed by our binaries. The excess grows to ~ 30 at i? prop ~ 10 h~ x kpc, 
and provides compelling evidence that the quasar autocorrelation function gets progressively steeper 
on sub-Mpc scales. This small scale excess can likely be attributed to dissipative interaction events 
which trigger quasar activity in rich environments. Recent small scale measurements of galaxy cluster- 
ing and quasar-galaxy clustering are reviewed and discussed in relation to our measurement of small 
scale quasar clustering. 

Subject headings: general - quasars: general - cosmology: general - surveys: observations - large-scale 
structure of the Universe 



1. INTRODUCTION 

A fundamental problem for cosmologists is to under- 
stand how quasars are embedded in the galaxy forma- 
tion hierarchy and to relate them to the gravitational 
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evolution of the structure of the underlying dark mat- 
ter. In the current paradigm, every massive galaxy is 
thought to have undergone a luminous quasar phase, and 
quasars at high redshift are the progenitors of the lo- 
cal dormant supermassive black hole population found in 
the cent ers of nearly all nearby bulge-dominated galax- 
ies ('e.g. iSmall fc Blandfordlll992UYu fc Tremainell2002ft 
. This fundamental connection is supported by the tight 
correlations between the masses of central black holes 
and t he velocity dispersions of their old stellar popula- 
tions llMagorrian et aT]ll998t iFerrarese fc Merritd l2000l 
iGebhardt et al.ll200Ct iTremaine et al.ll2002|) and bv com- 
paring the number density of black holes in the local Uni- 
verse to the l uminosity density pr oduce d by quasars at 
high redshift llSmall fc Blandfordlll992t lYu fc Tremainel 
2002). 

Qu asars are likely to reside in massive hosts ({Turned 
119911) and it has been suggested that they occupy the 
rarest peaks i n the initial Gauss i an den s ity fluctuation 
distrib ution llEfstathiou fc Reed 119881: iCole fc Kaisep 
-.891 iNusser fc Silkl Il993t IDiorgovski et all 119991 
Diorgovskil Il999t iDiorgovski et alJ 120031 iStiavelli et alJ 
2005). It is also thought that quasar activity is trig- 
gered by the frequent mergers which are a generic 
consequence o f bott o m up structure formation models 
(iBahcall et all Il997t iCarlberd 119901: lHaehnelt fc Reesl 
119931 iWvithe fc Loebll2002blh Both of these hypothesis 
imply that quasars should be highly biased tracers of the 
dark matter distribution: rare pea ks in the den sity field 
are intrinsically strongly clustered l|Kaiserlll984fl and the 
frequency of merger s is higher in dense environments 
ijLacev fc Coldir993j) . Measurements of quasar cluster- 
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ing can thus teach us about the environments of quasars 
as well as give clues to the dynamical processes which 
trigger quasar activity. Furthermore, a comparison of 
quasar clustering with the quasar luminosity function 
can be used to constrain the mean quasar lifetime 
l|Haiman fc Huil2001tlMart,ini fc Weinberdl2001fl as well 
as the relationship between the mass of central black 
holes and the circular velocities of their host dark halos 
<|Wvithe fc Loebll200^ . 

There have been many attempts to measure quasar 
cluste rin g, beginning w ith the pioneering work of lOsmeil 
(|1981l) . IShaverl l|1984D first detected quasar clustering 
using a clever technique to measure correlations from in- 
homogeneous catalogs and discovered that quasars were 
clustered similar to g alaxies in th e loca l universe, a re- 
sult later confirmed bvlKrusz ewskil (11 9881). Most recently, 
ICroom et all (|2001|) . |Porciani. Magliocchetti. fc Norberd 
(|2004^ (henceforth PMN). and ICroom et all <|2005|) meaT 
sured the clustering of quasars in the redshift range 
z = 0.3 — 2.2 from the ~ 15, 000 q uasars in the Two De - 
gree Field Quasar Survey (2QZ) <|Croom et al.ll2004b|) . 
They both find good agreement with a power law cor- 
relation function £(r) = (r/ro) -7 on scales r = 1 — 
35 hr 1 Mpc, with correlation length r$ ~ 5 hr 1 Mpc (co- 
moving) and slope 7 ~ 1.5, with only a weak dep endence 
on redshift and luminosity ijCroom et alJl200 2 j) . This 
agree s with previous measuremen t s lllovino fc S haver 
19881 lAndreani fc Cristiaml Il99l IMo fc Fand 119931 
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similar to the clustering of nearby galaxies. 

At redshift z > 2.5, quasars are rarer and thus quasar 
clustering has been much more difficult to measure. How- 
ever the mere existence of a few high redshift quasar 
pairs provides circumstantial evidence that quasars may 
have been much more highly clustered in the past. In 
the P alomar Transit Grism Survey of iSchneider et al.l 
( 1994), three quasar pairs with z > 3 and comoving sepa- 
rations 5 — 10 h~ x Mpc were found in a complete sample 
of 90 quasars covering 61.5 deg 2 . Anal ysis of cl uster- 
ing in this high redsh ift sample by IKundid (|1997|) and 
IStephens et al.l l|1997|) detected a statistically significant 
clustering signal, dominated by the three pairs, which 
implied a comoving correlation length tq ~ 50/i -1 Mpc. 
This is much larger than the correlation length of present 
day galaxies or z ~ 1.5 quasars. The only sub-arcminutc 
high redshift quasar pair known is a 33 " pair of quasars 
at z = 4.25 discovered serendipitously bv Sch neider et all 
(2000). Based on the discovery of this one object with 
proper transverse separation of 162 h~ x kpc, they es- 
timated the correlation length could be as large as ~ 
30 h~ x Mpc. iDiorgovski et alJ l)2003|) discovered a com- 
panion at z = 5.02 separated by 196" fro m the high 
redshi ft quasar at z = 4.96 discovered by iFan et all 
(1999), corresponding to a proper transverse separation 
of 896 hT 1 kpc. This is the highest redshift pair of 
quasars known. 

Even in the large quas ar sample studied by 
ICroom et all pfiol . PMN, and lCroom eTaf] pOfil . the 
smallest scale at which the correlation function can be 
measured is ~ 1 hr 1 Mpc. The reason for this is twofold. 
First, close quasar pairs with angular separations < 60", 
corresponding to ~ \h~ x Mpc at z ~ 1.5, are extremely 
rare, simply because at small separations, the correla- 
tion function does not increase as fast as the volume 



decreases. Second, because of the finite size of the op- 
tical fibers of the 2dF multi-object spectrograph, only 
one quasar in a close pair with separation < 30" can 
be observed. This limitation, referred to as a fiber colli- 
sion, is also a pro blem for the Sloan Digital Sky Survey 
llYork et al-lpQQO j SPSS) , for which this angular scale is 
55" JBlanton et alJl200^ . 

A significant motivation for studying small scale quasar 
clustering, is the existence of controversial population 
of quasar pairs discovered in the search for gravita- 
tiona l ly lensed quasars llKoc hanek. Falco. fc Munozl 
ll999HMortlock. Webster, fc Francis! 1 19991) . These close 
pairs have similar optical spectra, small velocity dif- 
ferences, and typically have separations in the range 
2" < A8 < 10" characteristic of group or cluster scale 
lenses. However, deep imaging shows no identifiable 
lenses in the foreground. Although a handful of wide 
separati on (A9 > 3") gravitational lenses hav e been dis- 
covered l|Walsh. CarswelL fc Wev mannlll979|). especially 
recen t ly in the SPSS ijlnada et alJ 12003 iQguri et alJ 
12004 l2005|) . the expected number of quasars lensed 
by groups and clusters is too small to account for 
all of the controver sial pairs (|Qguri fc Keetonl l200l 
iHennawi et alJ | 2005l). The poster child example is 
Q 2345+007 fWeedman et alJ 11982). the famous pair 
of z = 2.16 quasars with 7.1" separation. Although 
a plethora of papers can be found in the litcra- 
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of aga inst llPhinnev fc Blandford 1986; Piorgovski 
iSchneiderl Il993t IK ochanek. Falco. fc Munoz 
Mortlock. Webster. & Francis 1999i iCreen et al 



[19911 

2002) the lensing hypothesis for this system, the most 
compelling arg ument is based on the recent Chandra 
observations of ICreen et alJ l)2002fl who failed to detect 
diffuse X-ray emission associated with the potential lens. 

This population has thus led to much speculation 
about exotic mass concentrations which could be respon- 
sible for the apparent multiple imaging. It has been sug- 
gested that the lens es in these systems are 'dark' galaxies 
or galaxy clusters iSub ramaniam Rees. fc Chitra |1987t 
Duncarj l 99lllHawkins et al.lll997|lMalhotra et al.lll997 



: ^t"^iJ ; ^jpjJ^lJ|JJ| ; ^ ^11 ,A,±.MM^.,., II 

Peng et alll999tlKooDmans et alJl2000l) . that the y could 
be len sed by free floating ~ 10 14 M Q black holes l|Turneil 
Il991[) . or that they mig ht be instances of gravitational 
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lensing by cosmic strings l|Vilenki 
Hog an fc Naravanlll98l) . 

A much more plausible explanation is 
these controversial pairs are bin a ries rather 
lenses (Phinncv & Blandford 1986; Piorgovski 
Schneider 1993; Kochanek. Falco, fc Munozl 
Mortlock. Webster, fc Francis! [1999) and hence just 
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a manifestation of quasar clustering on small scales. 
iPiorgovskil 1)19911) first pointed out that this interpreta- 
tion implies a factor of ~ 100 more binary quasars over 
what is naively expected from extrapolating the quasar 
correlation function power law down to comoving scales 
< 100 h^ 1 kpc, and he proposed that this was due 
to the enhancement of quasar activity during merger 
events. Based on two clos e pairs found in the LBQS 
survey, iHewett et, al.l l)1998f ) similarly claimed an excess 
of ~ 100 over the expecta t ion fr om quasar clustering. 
IKochanek. Falco. fc Munozl l|1999D compared the optical 
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and radio properties of the controversial quasar pairs, 
and argued that they were all binary quasars, and 
similarly claimed that the excess binaries could be 
explained in a merger scenario. 

The study of binary quasars and small scale quasar 
clustering has been hindered by the small number of 
known examples and the heterogeneous mix of detection 
methods. In this paper we conduct a systematic search 
for binary quasars in the SDSS and 2QZ quasar sam- 
ples. We present a sample of 218 new binary quasars 
with proper transverse separations i? pr0 p < 1 h^ 1 Mpc, 
24 of which have angular separations < 10" correspond- 
ing to transverse proper separations i? pr0 p < 50/i _1 kpc, 
more than doubling the number of such systems known. 
A sub-sample of pairs selected with well defined crite- 
ria is constructed and we quantify its selection function. 
Based on this sample, we present the first measurement 
of the correlation function of quasars on the small scales 
10 h^ 1 kpc < i?p r0 p < 1 h~ x Mpc. We detect excess 
small scale clustering compared to the expectation from 
an extrapolation of the larger scale two point correlation 
function power law slope. 

The outline of this paper is as follows. In we dis- 
cuss color-selection criteria used to find binary quasars 
and describe the follow-up observations required to con- 
firm quasar pair candidates in § 0] Our binary quasar 
sample is presented in [JSJ We show that the number of 
binary quasars discovered thus far in the SDSS imply an 
excess of small scale quasar clustering in Sj^Jand we com- 
pare this result to small scale galaxy and quasar-galaxy 
clustering in HZj We summarize and conclude in SjSJ In 
the Appendix, we present tables summarizing the results 
of all of our follow-up observations, as well as a catalog 
of projected quasar pairs from the SDSS. 

Throughout this paper w e use the best fit W MAP 
(only) cosmological model of Spcrgc l et alJ l)2003f >. with 
n m = 0.270, n A = 0.73, h = 0.72. Because both 
proper and comoving distances are used, we will always 
indicate the former as i? pr0 p. It is helpful to remem- 
ber that in the chosen cosmology, for a typical quasar 
redshift of z = 1.5, an angular separation of Ad = 1" 
corresponds to a proper (comoving) transverse separa- 
tion of i? prop = 6 h~ x kpc (R — 15 h^ 1 kpc), and a 
velocity difference of 1000 km s _1 at this redshift cor- 
responds to a proper radial redshift space distance of 
Sprop = 4.4 h~ x Mpc (comoving s = 11 h~ x Mpc). 

2. QUASAR SAMPLES 

In this section we present a variety of techniques used 
to select quasar pair candidates. First, we describe the 
quasar catalogs which served as the parent samples for 
our quasar pair search. Then we introduce a statistic 
which quantifies the color similarity of two quasars. Fi- 
nally, we discuss each selection method in detail and de- 
scribe our follow-up observations. 

2.1. The SDSS Spectroscopic Quasar Sample 

The Sloan Digital Sky Survey uses a de dicated 2.5m 
telesc ope and a large format CCD camera ijGunn et al.l 
1998) at the Apache Point Observatory in New Mex- 
ico to obtain images in five broad bands (u, g, r, i and 
z, centere d at 3551, 4686, 6166 , 7480 and 8932 A, re- 
spectively; iFukugita et"aflll996t iStoughton et all 12002(1 



of high Galactic latitude sky in the Northern Galac- 
tic Cap. The ima ging data are pr ocessed by the as- 
trometric pipeline (iPier et all |2003|) and photometric 
pipeline ijLupton et all 12001)) . and a re photometrically 
calibrated to a st andard star network l)Smith et al.ll2002t 
iHogg et al.ll200lD . Additi onal details on th e SPSS data 
produ cts can be found in lAbazaiian et al.l l)2003l 120041 
2005). 

Based on this imaging data, spectroscopic targets cho- 
sen by various selection algorithms (i.e. quasars, galax- 
ies, stars, serendipity) are observed with two double 
spectrographs producing spectra covering 3800-9200 A 
with a spectral resolution ranging from 1800 to 2100. 
De tails of the s pectro sc opic observations c an be found 
in lYork et alJ fcOOOft iCastander et alJ 1)20011) . and 
IStoughton et all 1)20021) . A d iscussion of qu a sar tar - 
get selection can be found in iRichards" etatl l)2002a)) . 
The Third Dat a Release Quasa r Catalog contains 46,420 
quasars l)Schneider et all 120051) . Here, we use a larger 
sample of quasars, as we include non-public data: our 
parent sample includes 67,385 quasars with z > 0.3, of 
which 52,279 quasars lie in the redshift range 0.7 < z < 
3.0. Note also that we have used the Princeton/MIT 
spectroscopic reductions 17 which differ slightly from the 
official SDSS data release. 

Most quasar candidates are selected based on their lo- 
cation in multidimensional SDSS color-space. All mag- 
nitudes are reddening corrected following th e prescrip- 
tion in Schlegel. Finkb einer. fc Davis! (^98). Objects 
with colors that place them outside of the stellar lo- 
cus which do not inhabit specific "exclusion" regions 
(e.g., places dominated by white dwarfs, A stars, and M 
star-white dwarf pairs) are identified as primary quasar 
candidates. An i magnitude limit of 19.1 is imposed 
for candidates whose colors indicate a probable red- 
shift of less than rj 3; high-redshift candidates are ac- 
cepted if i < 20.2. Over 90% of SDSS-selected quasars 
follow a remarkably tight color-redshift r elation in the 
SDSS color-system l)Richards et alll2001a)) . In addition 
to the multicolor selection, unresolved objects brighter 
than i = 19.1 that lie within 1.5" of a FIRST radio source 
l)Becker. White, fc Helfandl Il995|) are also identified as 
primary quasar candidates. 

Supplementing the primary quasar sample described 
above are quasars targeted by other SDSS target selec- 
tion packages: Galaxy (the SPSS main and extended 
galax y samples lEisenstein et all 120011 iStrauss et all 
12002)) . X-ray (object s near the position of a ROSAT All- 
Sky Survey source. [Anderson et all 12003]) . Star (point 
source with unusual color), or Serendipity (unusual color 
or FIRST matches). No attempt at completeness is 
made for the last three categories; objects selected by 
these algorithms are observed if a given spectroscopic 
plate has fibers remaining after all of the high-priority 
classes (galaxies, quasars, and sky and spectrophotomet- 
ric calibrations) in t he field have been assigned fibers (see 
iBlanton et al.l l2003). Most of the quasars that fall below 
the magnitude limits of the quasar survey were selected 
by the serendipity algorithm. 

As we described in the introduction, the SDSS spec- 
troscopic survey selects against close pairs of quasars be- 
cause of fiber collisions. The finite size of optical fibers 

17 Available at http://spectro.princeton.edu 
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implies only one quasar in a pair with separation < 55" 
can receive a fiber. Follow-up spectroscopy is thus re- 
quired to discover quasar pairs. An exception to this 
rule exists for a fraction (~ 30%) of the area of the spec- 
troscopic survey covered by overlapping plates. For these 
plates the same area of sky was observed spectroscopi- 
cally on more than one occasion so that there is no fiber 
collision limitation. 

2.2. The SDSS Faint Photometric Quasar Sample 

iRichards et alJ (|2004) have demonstrated that faint 
(i < 21) photometric samples of quasars can be con- 
structed from the SDSS photometry, by separating 
quasars from stars using knowledge of their relative den- 
sities in color space. Each member of this catalog is 
assigned a probability of being a quasar, a photomet- 
ric redshift, and a probability that the p hotometric red- 
shift is correct (see IRichards et aT1l2004l for details). Wc 
searched for (and found) quasar pairs in a photomet- 
ric sample of 273,287 quasar candidates. Note that the 
faint photometric quasar used he re is based on the large r 
SDSS Data Release 3 area f DR3:IAbazaiian et al.l l20051. 
whereas that published in ( Richards ^^^JJ 2004) covers 
the smaller SDSS DR1 area l|Abazaiian et al J I2003) . 

2.3. The SDSS+2QZ Quasar Sample 

The 2dF Quasar Redshift Survey (2QZ) is a homo- 
geneous spectroscopic catalo g of 44,576 s tellar objects 
with 18.25 < bj < 20.85 flCroom et alJ 12004b) . Of 
these, 23,338 are quasars spanning the redshift range 
0-3 ^$ z ^ 2.9. Selection of quasar candidates is based on 
broad band colors (ubjr) from automated plate measure- 
ments of the United Kingdom Schmidt Telescope photo- 
graphic plates. Spectroscopic observations were carried 
out with the 2dF instrument, which is a multi-object 
spectrograph at the Anglo-Australian Telescope. The 
2QZ covers a total area of 721.6 deg2 arranged in two 
75° x 5° strips across the South Galactic Cap (SGP strip), 
centered on S = —30°, and North Galactic Cap (NGP 
strip, or equatorial strip), centered at 5 = 0°. The NGP 
overlaps the SDSS footprint, corresponding to roughly 
half of the 2QZ area. 

By combining the SDSS quasar catalog with 2QZ 
quasars in the NGP which have matching SDSS photom- 
etry, we arrive at a combined sample of 75,579 quasars 
with z > 0.3, of which 67,385 are from the SDSS and 8194 
from the 2QZ. For the clustering analysis in §Elwe will 
restrict attention to the redshift range 0.7 < z < 3.0, for 
which we define a combined SDSS/2QZ sample of 59,608 
quasars, of which 52,279 are from the SDSS and 6879 are 
from the 2QZ. 

3. QUASAR PAIR SELECTION 

Several different techniques are used to find binary 
quasars. For small separation pairs, A9 < 3", charac- 
teristic of the majority of gravitational lenses, binary 
quasars were discovered in the SD SS search for gravi- 
tationally lensed quasars (e.g., iQguri et al J 12005). For 
wider separations, both components are resolved and we 
can exploit the precise digital photometry of the SDSS 
to color select quasar pair candidates. Finally, quasar 
pairs can be found directly from the spectroscopy: the 
SDSS contains a fair number of overlapping spectroscopic 
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Fig. 1. — Distribution of the \ 2 statistic for the 64621 unique 
pair combinations of 359 quasars in the SDSS sample in the redshift 
interval 2.4 < z < 2.45. The median value of this distribution 
is 33.1, so that a quasar pair survey which aims to achieve 50% 
completeness in this redshift interval would have to observe all 
quasar pair candidates with \ 2 < 33.1. 



plates for which fiber collision does not limit the pair sep- 
aration and quasar pairs can be found by searching the 
SDSS+2QZ quasar catalog over regions where the survey 
areas overlap. 

3.1. Lens Selection 

For small separation pairs A8 < 3", the two images are 
unresolved or marginally resolved, as the SDSS imaging 
has a median seeing of 1.4". Candidates are selected by 
fitting a multi-component PSF model to a tlas images of 
each o f the SPSS quasa r s as d escribed in IPindor et alJ 
l)2003|) and llnada et alJ l)2003|) . We restricted atten- 
tion to candidates in the redshift range 0.7 < z < 3.0. 
Quasars with z < 0.7 are unlikely to be gravitational 
lenses, and the PSF fitting is complicated by the pres- 
ence of resolved host galaxy emission. The SDSS is bi- 
ased against gravitational lenses with z > 3.0 because 
the target selection algorithm for high-z quasars does not 
target objects classified as extended by the photometric 
pipeline, and most candidate le nses and binar i es app ear 
extended. See the discussion in IPindor et al.1 l)2003ft for 
more details. The number of quasars in this redshift 
range searched with our lens algorithm was 39,142, which 
makes up the parent sample of our lens search. This is a 
subset of the total number (52,279) of SDSS quasars in 
this range, as the lens algorithm was run on a sample of 
quasars defined at an earlier date. We refer to objects 
selected by this algorithm as the 'lens' sample. 

Follow up observations are required to confirm that the 
companion object is indeed a quasar at the same red- 
shift. The limiting magnitude for the companion objects 
is i < 21.0 (here and throughout we always quote red den- 
ing corrected asinh magnitudes l|Lupton et al]ll999|) '). as 
fainter objects are too difficult to observe from the 3.5m 
telescope at Apache Point Observatory (APO), where 
most of the follow-up observations were conducted, even 
in the best conditions (see §0J. 

3.2. x 2 
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Fig. 2. — Completeness of x 2 statistic as a function of redshift for 
a survey which observes all quasar pair candidates with x 2 < 20. 
Dispersion in the color-redshift relation of quasars ( Richards et a),. 
I2001al) gives rise to a broad distribution of \ 2 for pairs of quasars 
at the same redshift (see Figure 0. Restricting the observations of 
pair candidates to those with x 2 < 20 results in ~ 35% complete- 
ness in the redshift range 0.70 < z < 3.0. 



Although quasars have a wide range of luminosities, 
the majority have s imilar optical /ultravio let spectral en- 
ergy distributions. ( Ri chards et al.ll2001a|) demonstrated 
that most quasars follow a relatively tight color-redshift 
relation in the SDSS filter system; a property which 
has been e xploited to calculate photometric redshifts 
of quasars (Ric hards et alJl2001bt iBudavari et alJl200H 
IWeinstein et al.N2004j) . It is thus possible to efficiently 
select pairs of quasars at the same redshift by searching 
for pairs of objects with similar, quasar-like colors. 

To this end, we define a statistic that quantifies the 
likelihood that two astronomical objects have the same 
color. Recall that a color u — g, is a statement about 
flux ratios, f u /f 9 . Thus if two objects have the same 
color, then their fluxes should be proportional. Given 
the fluxes /{" of the first, we can ask whether the fluxes 
of the second are consistent with the proportionality 



f? = a /r 



where f m is a five dimensional vector of fluxes (one for 
each SDSS band), m designates the filter, and this rela- 
tionship holds with the same proportionality constant A 
in all bands. 

The maximum likelihood value of the parameter A, 
given the fluxes of both objects, can be determined by 
minimizing the \ 2 
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(2) 



where we have dropped a term corresponding to the nor- 
malization of the likelihood because of its slow variation 
with the parameter A. The er m are the photometric mea- 
surement errors, but do not include the intrinsic scatter 
about the mean color-redshift relation (see discussion be- 
low). We thus arrive at the implicit equation for A 



which can be solved in a few iterations. This value is 
inserted into eqn. Q), reducing the number of degrees of 
freedom in the \ 2 to four, i.e. the number of independent 
colors one could have formed from the five magnitudes. 

If the fluxes of the two objects are proportional and 
if the photometric errors are distributed normally, this 
statistic follows the chi-squared distribution with four 
degrees of freedom, and the typical value will be \ 2 ~ 4. 
However the colors of two quasars at the same redshift, 
although similar, will in general not be exactly propor- 
tional. Fluctuatio ns about the median color-redshift re- 
lation of quasars {Richards et al]l2f)01a() will result in an 
additional source of 'dispersion' in our color similarity 
statistic. This leads to a much broader distribution of 
X 2 than expected from Gaussian statistical errors. Fig- 
ure shows the distribution of \ 2 f° r the 64621 unique 
pair combinations of 359 quasars in in the redshift inter- 
val 2.4 < z < 2.45. The median value of this distribution 
is 33.1, so that a quasar pair survey which aims to achieve 
50% completeness in this redshift interval would have to 
observe all quasar pair candidates with \ 2 < 33.1. Also 
notice that a long tail in this distribution extends even 
beyond \ 2 ~ 100 because of ou tliers from the media n 
color-redshift relation of quasars ({Richards et alJ 2001a1. 
caused by broad absorption line features or reddening 
{Richards et al.ll2003UHoDkins et alJl2004jl . 

Our survey for close pairs of quasars thus involves a 
tradeoff between completeness and efficiency, since toler- 
ating larger values of % 2 will increase the number of false 
pair candidates. Our follow up observations targeted 
quasar pair candidates with x 2 < 20. This threshold im- 
plies a certain level of completeness for our survey, which 
varies redshift, as the dispersion in th e color-redshift 
relatio n of quasars depends on redshift (Ric hards et alJ 
l2001a|L We quantify this incompleteness by dividing the 
SDSS quasar sample into redshift bins of Sz — 0.043, 
and computing the fraction of all the unique combina- 
tions of pairs in each bin with x 2 < 20. Figure [21 shows 
the completeness of our quasar pair survey as a function 
of redshift. For the redshift range 0.70 < z < 3.0, where 
most of our binary quasars lie, the \ 2 < 20 cut results in 
~ 35% completeness. 

Finally, note that the statistic defined by eqn. J5J) uses 
an isotropic 'metric' in color-space. This would not be 
the case had we included the variance of the color-redshift 
relation a m (z) in our errors o - ^„,, similar to the pro- 
cedure used by iRichards et alJ pOOlblk IWeinstein et 
( 2004) to determine photometric redshifts of quasars. In 
retrospect, this would be a more suitable procedure for 
finding binary quasars. However, we were also conduct- 
ing a search for wide separation gravitational lenses, and 
for lenses there is no color-redshift scatter (after all, it is 
the same quasar observed twice!). 

We applied our color similarity statistic to all objects 
within the annulus 3" < A9 < 60" around the 59,608 
quasars (0.7 < z < 3.0) in the combined SDSS+2QZ 
catalog. We refer to the binaries selected by searching 
around the SDSS+2QZ quasars as members of our 'x 2 ' 
sample. To be considered for follow up observations, a 
quasar pair had to meet the following criteria 
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i < 21.0 (4) 
a <0.2. 

The minimum redshift was imposed because of our 
desire to find wide separation g ravitational lenses 
ijlnada et alJl2003t lOeuri et alJ l2004). The a l < 0.2 re- 
quirement gets rid of objects with very large photometric 
errors due to problems with dcblcnding or very poor im- 
age quality. In addition, we required that the companion 
objects be optically unresolved in the SDSS imaging to 
avoid contamination from galaxies. Nearly all quasars at 
z > 0.70 should be unresolved in the SDSS imaging, with 
the exception of very small separation pairs < 3", how- 
ever, these pair candidates are selected by our algorithm 
described in § 13.11 

We also used the criteria in eqn. JSJ to search for pair 
candidates in the catalog of 273,287 faint photometric 
quasars. We restricted attention to members of the cat- 
alog only, and did not consider other nearby photometric 
objects. We will refer to binaries discovered in this cat- 
alog as our 'photometric sample'. The same criteria as 
in eqn © were used, but there were no lower or upper 
limits on redshift since spectroscopic redshifts are not 
available for this sample. 

3.3. Overlap and Spectroscopic Selection 

As mentioned previously, the SDSS contains a fair 
number of overlapping spectroscopic plates for which 
fiber collision does not limit the pair separation. Fur- 
thermore, quasar pairs can be found below the fiber col- 
lision limits by searching the SDSS+2QZ quasar catalog 
over regions where the survey areas overlap. Finally for 
separations 9 > 60", larger than the SDSS fiber collision 
scale 18 , quasar pairs can be found in the entire SDSS 
area. We refer to pairs found in the spectroscopic cat- 
alog with 9 < 60" as our 'overlap' sample and we refer 
to those with 9 > 60" as our 'spectroscopic' sample. No 
color similarity criteria were applied to these objects, and 
the only magnitude or err or limits are those imp osed by 
the SDSS target selection {Richards et alJl2002aj) . 

3.4. Previously Known Binaries 

It is instructive to ask whether previously known bi- 
nary quasars are selected by our selection techniques. In 
Table^we list the nine previously known binary quasars 
with 0.7 < z < 3.0 which are in the SDSS footprint, of 
which seven were recovered. The column labeled 'sample' 
indicates the sample for which each binary was selected 
as a candidate. Six of these binaries were listed in the 
compilation of binary quasars on the CASTLES 19 web- 
site. The others are SDSS J2336-0107 fGregg et al.l2002(l 
and SDSS J1120+6711 (|Pindor et alJl2005j) . discovered 
recently in the SDSS search for gravitational lenses, and 
2QZ J1435 +0008, the 33" pair of quasars discovered in 
the 2QZ bv [Miller et all p004|) . 

Of the nine binaries listed in the table, all but 
Q 1343+2650, Q 1635+267, and 2QZ J1435+0008 had 
at least one member of the pair in the SDSS spectro- 
scopic sample. Neither member of the pairs Q 1343+2650 

18 Although the fiber collision limit is 55", for operational pur- 
poses we take it to be 60" to give a small buffer from the actual 
limit where the fiber tiling may still be imperfect. 

19 Available at http://cfa-wwww.harvard.edu/castles 



and 2QZ 1435+0008 were targeted for spectroscopy be- 
cause these quasars are below the flux limit of the SDSS 
quasar catalog (i < 19.1 for quasars in this region of color 
space). The brighter member of Q 1635+267 is above 
the flux limit, but this area of sky has only been imaged 
and has yet to be spectroscopically observed. However, 
both members of this pair are members of the faint pho- 
tometric catalog, and indeed, this pair was selected as 
a candidate by our photometric selection. Both m em- 
bers of the binary HS1216+5032 ijHagen et al.lll996|) re- 
ceived SDSS fibers, so that this binary is a member of 
our overlap sample. The brighter of the two members 
of the famous double quasar Q 2345+0007 received an 
SDSS fiber, and this pair was selected as a pair candi- 
date by our \ 2 algorithm. Of the two binaries which 
were not recovered, Q 1343+2650 was missed because 
it was below t he SDSS flux limits and the quasar pair 
Q 1120+0195 ijMevlan k Diorgovskil fl989^ was missed 
because its x 2 = 20.6, is just above the cutoff \ 2 < 20 of 
the 'x 2 ' selection algorithm. 

4. SPECTROSCOPIC OBSERVATIONS 

Candidates in our lens, x 2 , and photometric samples 
require follow up spectroscopy to confirm the quasar pair 
hypothesis, which we describe in this section. 

The SDSS images of the candidate quasar pairs and the 
spectrum of the quasar with an SDSS or 2QZ spectrum 
were visually inspected to reject bad imaging data and 
possible spectroscopic misidentifications. Color-color di- 
agrams for each candidate were also visually inspected, 
and pairs for which the companion object overlaps the 
stellar locus (see e.g. iRichards et al.ll2001aj) were given a 
lower priority. 

The result of a successful follow up observation of a 
quasar pair candidate falls into one of four categories: (1) 
a quasar-quasar pair at the same redshift (2) a projected 
pair of quasars at different redshifts (3) a quasar-star 
pair (4) a star-star pair (for the photometric catalog). 
As an operational definition, we consider quasar pairs 
with velocity differences of \Av\ < 2000 km s -1 to be at 
the same redshift, since this brackets the range of veloc- 
ity differences caused by both peculiar velocities, which 
could be as large as ~ 500 km s -1 if binary quasars 
reside in rich environments, and redshift uncertainties 
caused by blueshifted broad lines (~ 1500 km s -1 ) 
{Richards et al.ll2002b|h 

Spectra of the vast majority of our quasar pair can- 
didates were obtained with the Astrophysical Research 
Consortium (ARC) 3.5m telescope at the Apache Point 
Observatory (APO), during a number of nights between 
March 2003 and January 2005. In addition, two of the 
binary quasars in our sample were confirmed at other 
telescopes: SDSSJ1600+0000 was confirmed at the ESO 
3.58m New Technology Telescope, and SDSSJ1028+3929 
was discovered at the Hobby-Eberly Telescope. Higher 
signal to noise ratio spectra of five of the binary quasars 
in our sample were obtained at the 10m Keck I & Keck 
II telescopes. 

For the ARC 3.5m observations, we used the Double 
Imaging Spectrograph (DIS), a double spectrograph with 
a transition wavelength of 5350 A between the blue and 
red side. The observations were taken with low resolution 
gratings, with a dispersion of 2.4 A pixel -1 in the blue 
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TABLE 1 

Previously Known Binary Quasars in the SDSS Footprint 



Name 


z 


A9 


^prop 


|Ad| 


i\ 


ii 


x 2 


Sample 


Notes 


SDSSJ 1120+6711 


1.49 


1.5 


9.0 


<100 


18.49 


19.57 




lens 


1 


Q 1120+0195 


1.47 


6.5 


39.6 


630 


15.61 


20.26 


20.6 




2 


HS 1216+5032 


1.46 


9.1 


55.1 


50 


16.76 


18.31 


534.1 


overlap 


3 


Q 1343+2650 


2.03 


9.2 


55.4 


200 


19.14 


19.97 


113.7 




4 


LBQS 1429-008 


2.08 


5.1 


30.8 


260 


17.40 


20.74 


4.3 


x 2 


5 


2QZ 1435+0008 


2.38 


33.2 


194.8 


760 


19.90 


20.56 


7.0 


x 2 


6 


Q 1635+267 


1.96 


3.9 


23.2 


30 


19.04 


20.26 


7.7 


photo 


7 


SDSS J2336-0107 


1.29 


1.7 


10.0 


240 


19.26 


18.94 




lens 


8 


Q 2345+007 


2.16 


7.1 


42.4 


480 


18.68 


20.45 


0.3 


x 2 


9 



Note. — The redshift of the binary quasar is z, A8 is the angular separation, R is the 
proper transverse separation, |Ati| is the velocity difference between the two quasars in 
km s , i\ and %i are extinction corrected i-band magnitudes of the brighter and fainter 
quasar respectively, and \ 2 i s the value of our color similarity statistic, computed only for 
pairs with A9 > 3". The last column labeled 'Sample' indicates which of our selection 
al gorithms recovered the binary. 

1 (Pindor ct al. 2005) SDSS binary quasar 

2 IMevlan fc Diorgovskilll98flh Q 1120+019.-, is also named PHL 1222 and I'M lii 

3 (Hagen ct al. 1996). Spectra of both quasars are in the SDSS sample 

4 ICrampton et a l. 1988) Both quasars are below the flux limit of the SDSS quasar survey. 

5 IHewett et aLlll98Hh 

6 IMiller et al.1120041) 

7 ISramek fc Wccdmanl IT978T) This region of sky has been imaged by the SDSS but not 
yet spectroscopically observed. Both quasars are members of the faint photometric quasar 
sa mple 

8 IGregg et a3[2 002) SDSS binary quasar 

9 IWeedman et al.|[l982T) 



side and 2.3 A pixel -1 on the red side, and a resolution 
of roughly 2 pixels. A 1.5" slit was used and we oriented 
the slit at the position angle between the two quasars so 
that both members of the pair could be observed simulta- 
neously. The final spectrum covers the wavelength range 
of 3800 A to 10,000 A. The wavelength scale is calibrated 
with a polynomial fit to lines from an Ar-He-Ne lamp; 
the typical rms error of the fits is smaller than 0.5 A. 
Observatio ns of a variety of o ptical spectrophotometric 
standards (O ke & Oundll983l) provided flux calibration; 
however, most of the candidates were not observed un- 
der photometric conditions, nor were they observed with 
the slit oriented at the parallactic angle. Exposure times 
ranged from 1200 seconds for candidates with i ~ 18 to 
2400 seconds for the faintest candidates, i ~ 20.8. 

The binary quasar SDSSJ1600+0000 was discovered at 
the ESO New Technology Telescope 3.58m on UT 2001 
April 18, using the red CCD of the ESO Multi-Mode In- 
strument (EMMI). The #13 grating (2.66 A/pixel, R ~ 
600 at 6000 A) was used with an OG530 blue-blocking 
filter. Relative spectrophotometric cali bration wa s ob- 
tained through observations of GD 108 l)Okelll990(K but 
this calibration is uncertain blueward of 5350 A and red- 
ward of 9300 A. 

The binary SDSSJ1028+3939 was identified from ob- 
servations obtaine d using the M arcario Low Resolu- 
tion Spectrograph ( Hill et al.ll998l) on the Hobby- Eberly 
Telescope (HET) on UT 2005 January 14. The HET ob- 
servations were obtained with a 300 line mm -1 grating 
and GG385 blocking filter. Spectra of both components 
were obtained simultaneously using a 1.5" slit. The spec- 
tra covered the rang e 4400-8000 A at a resolution of 18 A. 



Higher signal to noise ratio spectra of both components 
of the quasar pairs SDSSJ0955+6045, SDSSJ1010+0416, 
and SDSSJ1719+2549 were obtained on UT 2003 Febru- 
ary 5-6, using the Ec helle Spectrograph and Imager (ESI; 
IEpps fc Millerll998|) on the Keck II telescope. The seeing 
was 0'.'6. The ESI has a dispersion of 0.15-0.3 A pixel" 1 
over a wavelength range of 4000 — 10500 A, and the 1" 
slit used for these observations projects to 6.5 pixels. The 
900s exposures were obtained with the slit aligned at the 
position angle of the components of the pair. 

Spectra of both quasars in the pair SDSSJ0248+0009 
were obtained on UT 1999 October 17 using the Low- 
Resolu tion Imaging S pectrograph Keck II telescope 
fT;RTS: l(Tke~et al.l1 995ft . which was before LRIS was com- 
missioned as a double spectrograph. The 300 line mm -1 
grating blazed at 5000 A was used, giving a spectral cov- 
erage of 5220 A and a dispersion of 2.55 A pixel" 1 . The 
0".7 longslit was used and the seeing was 0".6. A single 
900s exposure was obtained with the slit aligned at the 
position angle of the components of the pair. 

We obtained spectra of both quasars in the pair 
SDSSJ0048-1051 on UT 2003 September 27 using the 
LRIS Double Spectrograph on the Keck I telescope. The 
D560 dichroic was used, which splits light between the 
blue arm and the red arm at 5600 A. On the blue side, 
the 600 line mm -1 grism blazed at 4000 A was used, giv- 
ing a spectral coverage of 2590 A and a dispersion of 0.63 
A pixel" 1 . On the red side, the 400 line mm -1 grating 
was used blazed at 8500 A, giving a spectral coverage 
of 3810 A and a dispersion of 1.86 A pixel" 1 . The 1".0 
longslit was used and the seeing was 0."6. A single 900s 
exposure was obtained at an airmass of 1.2 with the slit 
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aligned at the position angle of the components of the 
pair. 

All the data were reduced using standard procedures in 
the IRAF 20 package, supplemented by IDL routines bor- 
rowed from the SDSS spectroscopic reduction software 
and adapted to the different instruments. Quasar red- 
shifts were determined by cross correlating the quasar 
spectra with the first four eigenspectra of a principal 
component decompo sition of the SDSS quasar sample 
llSchlegel et al.ll20f)l . 

5. BINARY QUASAR SAMPLE 

In this section we present a sample of 218 new quasar 
pairs with transverse separations -R pr0 p < 1 h Mpc 
over the redshift range 0.5 < z < 3.0. Of these, 65 
have angular separations 9 < 60", i.e. below the SDSS 
fiber collision scale. Our 26 new binaries with transverse 
proper separations i? pr0 p < 50 h^ 1 kpc (9 < 10") more 
than doubles the number of known binary quasars with 
separations this small. Table lists relevant quantities 
for 33 binaries with 3" < 9 < 60" discovered from our 
X 2 and photometric samples. The last column indicates 
which algorithm which was used to select the binary. The 
binaries with 9 < 3" discovered from our lens sample 
are shown in Table |3 and the overlap and spectroscopic 
binaries found by searching for pairs in the SDSS+2QZ 
quasar catalog are shown in Table 0] Table 03 gives a 
summary of the number of binary quasars selected by 
each algorithm described in § [21 For completeness, the 
Appendix includes tables of projected pairs of quasars at 
different redshifts, as well as projected quasar-star pairs. 

The distribution of redshifts and proper transverse sep- 
arations probed by these binary quasars is illustrated 
by the scatter plot in Figure The (magenta) up- 
side down triangle are members of the lens sample, the 
(green) squares are members of the photometric sample, 
and (red) triangle are in the x 2 sample. The (blue) open 
circles are members of the overlap sample and smaller 
(blue) dots are pairs in the spectroscopic sample. The 
dashed curve shows the proper transverse distance cor- 
responding to 9 = 3", which divides the lens sample from 
the other samples, and the dotted line indicates the dis- 
tance corresponding to 6 = 60", above the fiber colli- 
sion limit so that pairs can be found in the spectroscopic 
quasar catalog. It should be noted that the distribution 
of points in Figure [3] reflects some biases in our observa- 
tional program. In particular, we tended to observe small 
separation pairs first, and we were much more likely to 
observe candidates with z > 2, because quasar pairs at 
these redshifts are of i nterest for studying the Lya forest 
ijHennawi. et alJl2"005j) . 

We next discuss the possibility that some of the quasar 
pairs in this sample are strong gravitational lenses rather 
than binaries. After showing spectra of some of the more 
notable binaries, we define a sub-sample of binaries se- 
lected homogeneously which we will use in our analysis 
of small scale quasar clustering in § 

5.1. Contamination by Gravitational Lenses 

20 IRAF is distributed by the National Optical Astronomy Ob- 
servatories, which are operated by the Association of Universities 
for Research in Astronomy, Inc., under cooperative agreement with 
the National Science Foundation. 
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Fig. 3. — Range of redshifts and proper transverse separations 
probed by the binary quasars published in this work (see Tables 1^1 
m m and |5J . The blue circles are binary quasars identified in the 
SDSS spectroscopic sample. The region to the left of the dotted 
curve is excluded because of fiber collisions (8 < 60"), with the ex- 
ception of the binaries discovered from overlapping plates, which 
are indicated by the larger open blue circles. The magenta circles 
are members of the lens sample, red circles are from the x 2 sample, 
and green squares are binaries from the photometric sample. Be- 
cause of the fiber collisions, the vast majority of small separation 
pairs R < 300 h' 1 kpc were discovered from our follow-up ob- 
servations. The dashed curve indicates the transverse separation 
corresponding to 9 = 3" below which binaries are found with our 
lens algorithm. Although we publish only pairs with separations 
fiprop < 1 h~ 1 Mpc in this work, Figureyjshows all pairs in the 
SDSS+2QZ catalog out to 3 Mpc for the sake of illustration. 



It is possible that some of the quasar pairs with 
image splittings < 15" in our sample could be wide 
separation strong gravitational lenses rather than bi- 
nary quasars. Indeed, the recently discove red quadru- 
ply imaged lensed quasar SDSSJ1004+4112 <|Inada et alJ 
2003; Oe uri et alJ I2004D with a maximum image sep- 
aration of 14". 6 was discovered as part of our follow 
up campaign to discover quasar pairs, as w ere two new 
gravi tational lenses with separations > 3" ijQguri et alJ 
2005). We thus review the set of objective criteria 
that must be satisfie d for a quasar pair to be classi- 
fied as a binary or lens |K ochane k. Falco. &: Munodl999t 
Mortlock. We bster, fc Francis! 119991) and briefly discuss 
why we have concluded that the binary hypothesis is cor- 
rect for our pairs. 

A quasar pair can be positively confirmed as a bi- 
nary if th e spectra of the images are vastly differ- 
ent (c.f. ICxregg et all 12001 . if only one of the im- 
ages is radio-loud (an Q 2 R pai r, in the notation of 
iKochanek. Falco. fc Munod 119991) . or if the quasars 
hosts are detected and they are not clearly lensed. The 
sufficient conditions for a pair to be identified as a lens 
are the presence of more than two images in a lensing 
configuration, the measurement of a time delay between 
images, the detection of a plausible deflector, or the de- 
tection of lensed host galaxy emission. 

For the majority of pairs in our sample, the APO dis- 
covery spectra have too low a signal to noise ratio to make 
convincing arguments for or against the lensing hypoth- 
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esis based on spectral dissimilarity. The exceptions are 
the five binaries for which we have high signal to noise 
ratio Keck spectra (three from ESI and two from LRIS). 
We comment on the spectral similarity of these binaries 
below. 

Although the absence of a deflector in images of a 
quasar pair does not strictly speaking confirm the bi- 
nary hypothesis, it certainly makes it more plausible. 
For small separation (AO < 3") gravitational lenses, the 
lens galaxies are rarely detected in the relatively shallow 
SDSS imaging. However, for the wider separation lenses 
with AO > 3", like SDSSJ1004+4112 (A6» max = 14.62", 
z, BT1s = 0.68) or Q 0957+561 (A0 = 6.2", z lcns = 0.36; 
I Walsh. Carswell. fc Wevmannlll979|) . where the lens is a 
bright galaxy in a cluster or group, the lens galaxies are 
detected in the SDSS imaging, though it is quite possi- 
ble that fainter high redshift lens galaxies or clusters in 
other wide separation lens systems would go undetected. 

Using the University of Hawaii 2.2m telescope, we 
have taken deeper optical (i < 24; z < 22) or near 
infrared images (J < 22; H < 22) of all the binaries 
in our sample with separations < 4" and of a sub- 
set of those with wider separations. None of the images 
showed lens galaxies in the foreground. For the wider 
> 4" pairs, the lensing hypothesis would require a very 
bright massive galaxy in a group or cluster. In partic- 
ular, because the typical redshift range of our binaries 
is z = 1.5 — 2, the most probable lens redshifts would 
be in the range z — 0.3 — 0.7, which would likely have 
been detected in the SDSS imaging. Furthermore, these 
wide separation multiply i maged quasars are extremely 
rare ijHennawi et al]l2005|) . thus we are confident that 
the pairs in our sample are all binaries. 

5.2. Sample Spectra of Binaries 

Keck ESI spectra of the three binaries 
SDSSJ0955+6045, SDSSJ1010+0416, and 

SDSSJ1719+2549 are shown in Figure H The sig- 
nal to noise ratios of these spectra are high enough that 
we can make arguments against the lensing hypothesis 
based on spectral dissimilarity. In particular, for 
SDSSJ0955+6045 the narrow [OIII] emission lines are 
significantly stronger in one of the quasars than the 
other. In SDSSJ1010+0416, the CIII] emission lines 
have a velocity offset of ~ 2000 km s _1 , although 
this offset is less apparent in Mgll, which tends to 
be a better tracer of t he systemic redshifts of quasars 
lIRichards et al.ll2002bft . Finally, for SDSSJ1719+2549, 
the peak to continuum flux ratios of the Mgll broad line 
differ by a factor of ~ 1.5 between the two quasars. 

Figure shows Keck LRIS spectra of SDSSJ0048-1051 
(moderate resolution) and SDSSJ0248+0009 (low res- 
olution). For SDSSJ0048-1051 the profiles of all the 
emission lines differ significantly, especially Mgll. For 
SDSSJ0248+0009, the peak to continuum flux ratios dif- 
fer significantly for CIV, CIII], and Mgll. 

Figures El and [7| show SDSS and APO spectra of six 
other binaries in our sample. 

5.3. Clustering Sub-sample 

In this section, we define a statistical sub-sample of 
binary quasars which we will use to measure the quasar 
correlation function in In the various samples 



used to identify the binary quasars in Tables |21 01 and 
01 were described. These various selection algorithms se- 
lected quasar pairs over different angular scales, with dif- 
ferent limiting magnitudes, and varying degrees of com- 
pleteness. Here we combine these samples in a coher- 
ent way, which will allow us to quantify their selection 
function. We pay special attention to the completeness 
of each sample used and the parent sample of quasars 
searched to define each sample. 

For the clustering analysis we restrict attention to 
quasars in the redshift range 0.7 < z < 3.0 with ve- 
locity differences \Av\ < 2000 km s . We use the lens, 
X 2 , overlap, and spectroscopic samples. Our approach is 
to stitch the samples together as a function of angle. The 
photometric quasar pairs are not in the SDSS+2QZ cat- 
alog so the selection function and completeness of these 
binaries is more difficult to quantify. 

All pairs with < 3" come from the lens sample. The 
parent sample of quasars for this angular range is the 
39,142 SDSS quasars to which we applied the lens algo- 
rithm. The completeness is the pro duct of the intrinsic 
completeness of t he lens algorithm fPi ndor et al.ll2003l 
llnada et aUl2003h and the fraction of candidates which 
have had spectroscopic confirmation thus far in the sur- 
vey. 

For pairs in the range 3" < < 60" we use the x 2 sam- 
ple. Quasar pairs from the overlap sample which also 
meet the x 2 selection criteria in eqn. (0 are included, 
and can be thought of as follow up observations which 
came for 'free' from the overlapping plates. Of the 21 bi- 
naries with 0.7 < z < 3.0 in the overlap sample, 8 satisfy 
X 2 < 20 (see Table 0J, and are included in the cluster- 
ing sub-sample. The completeness of binary quasars in 
the range 3" < < 60" is the product of completeness 
of the x 2 selection, C(z|x 2 < 20), and the fraction of 
candidates observed thus far. The parent sample around 
which we searched with the x 2 algorithm is the combined 
SDSS+2QZ quasar sample of 59,608 quasars in the range 
0.7<z<3.0. 

For separations > 60" , we use pairs found in the 
SDSS spectroscopic catalog of 52,279 quasars. We re- 
strict attention to the SDSS (rather than SDSS+2QZ), 
because the completeness for detecting quasar compan- 
ions is very high if we restrict attention to compan- 
ions above the SDSS flux limit for low redshift quasars. 
iVanden Berk et all l)2005l) measured a completeness of 
~ 95% for quasars in the range 0.3 < z < 3.0 with 
i < 19.1. Thus we only include quasar pairs > 60" in 
our clustering sample provided that at least one member 
of the pair is brighter than this flux limit. 

Finally, any of the previously known binaries listed in 
Table^which satisfied any of the criteria for the cluster- 
ing sub-sample are also included. Thus we include the 
binaries SDSS J1120+6711 and SDSS J2336-0107 as part 
of the lens sample, and LBQS 1429-0008, Q 2345+007, 
and 2QZ 1435+0008 are included as part of our x 2 sam- 
ple. 

Of the 65 quasar pairs with angular separations < 
60" which we publish in this work, 35 are included in our 
clustering sub-sample along with five previously known 
binaries for a total of 40 sub-arcminute pairs. The dis- 
tribution of redshifts and proper transverse separations 
of our clustering sample is is illustrated by the scatter 



10 



Honnawi et al. 



TABLE 2 

Binary Quasars with separations 3" < AO < 60" Discovered From Follow-Up Observations 



Name 


RA (2000) 


Dec (2000) 


u 


< 


/ 


j 




i 


z 


SDSSJ0048-1051A 


00 


:48 


:00 


.77 


-10: 


:51 


:48 


.6 


20 


.99 


20 


.52 


20. 


.18 


19. 


.94 


19. 


.91 


APOJ0048-1051B 


00 


:48 


;00 


.96 


-10: 


:51 


:46. 


.2 


20 


.39 


20 


.04 


19 


.70 


19 


30 


19. 


.28 


SDSSJ0054-0946A 


00 


:54 


:08 


.47 


-09: 


:46: 


:38 


3 


18 


.17 


17 


.90 


17 


.71 


17 


.51 


17. 


.31 


APOJ0054-0946B 


00 


:54 


;08 


.04 


-09: 


16 


:25 


,7 


20 


.87 


20 


.70 


20 


37 


20. 


11 


19. 


,74 


SDSSJ0201+0032A 


02 


:01 


:43 


.49 


+00: 


:32: 


:22. 


,7 


19 


.99 


19 


.39 


19. 


47 


19. 


.41 


19. 


.19 


APOJ0201+0032B 


02 


:01: 


:42 


.25 


+00: 


:32 


:18. 


.5 


20 


.80 


20 


.30 


20 


.14 


20. 


.12 


19. 


92 


SDSSJ0248+0009A 


02 


:48 


:20 


.80 


+00: 


09 


56 


.7 


19 


.44 


19 


.24 


19 


.23 


18. 


98 


19. 


00 


APOJ0248+0009B 


02 


:48 


:21 


.26 


+00: 


09 


:57. 


3 


20 


.77 


20 


.71 


20. 


.74 


20 


57 


20. 


39 


APOJ0332-0722A 


03 


:32 


:38 


.38 


-07: 


:22 


.15 


.9 


20 


.24 


20 


.29 


20 


.00 


19 


.78 


19 


.63 


APOJ0332-0722B 


03 


:32 


:37 


.19 


-07: 


:22: 


19 


,6 


20 


.63 


20 


.57 


20 


.22 


19. 


.96 


19 


.70 


SDSSJ0846+2749A 


08 


:46 


:31 


.77 


+27: 


:49: 


:21. 


9 


19 


.55 


19 


.57 


19. 


66 


19 


.49 


19. 


.26 


APOJ0846+2749B 


08 


:46 


:30 


.38 


+27: 


:49 


:18. 


1 


19 


.82 


19 


.88 


19 


.82 


19. 


.71 


19. 


.55 


SDSSJ0939+5953A 


09 


:39 


:48 


.78 


+59: 


:53: 


48. 


,7 


20. 


.40 


19 


.85 


19. 


80 


19 


.79 


19 


.45 


APOJ0939+5953B 


09 


:39 


:46 


.56 


+59: 


53 


:20 


,7 


19 


.31 


18 


.67 


18 


.57 


18. 


56 


18. 


42 


SDSSJ0955+6045A 


09 


:55 


:24 


.37 


+60: 


:45: 


:51. 





20 


.84 


20 


.37 


20. 


29 


20 


.25 


20 


.34 


APOJ0955+6045B 


09 


:55 


:25 


.37 


+60: 


:45 


33. 


.8 


20 


.68 


20 


.65 


20 


62 


20 


.72 


20 


30 


APOJ0959+5449A 


09 


:59 


:07 


.46 


+54: 


:49: 


:06 


.7 


20 


.26 


20 


.06 


19. 


95 


19 


.76 


19 


.73 


APOJ0959+5449B 


09 


:59 


:07 


.05 


+54: 


:49: 


:08. 


4 


20. 


.48 


20 


.61 


20. 


.43 


20 


.28 


19 


.82 


SDSSJ1010+0416A 


10 


:10 


:04 


.98 


+04: 


:16: 


36 


,2 


20 


.18 


20 


.08 


20. 


.00 


19 


.87 


20. 


.08 


APOJ1010+0416B 


10 


:10 


:04 


.37 


+04: 


:16: 


:21. 


6 


20 


.29 


20 


.19 


20. 


03 


19 


.87 


19. 


,94 


SDSSJ1014+0920A 


10 


:14 


:11 


.43 


+09: 


:20: 


17. 


,7 


19 


.95 


19 


.14 


19. 


12 


18 


.92 


18 


.73 


APOJ1014+0921B 


10 


:14 


;10 


.29 


+09: 


21 


:01. 


.7 


20 


.81 


20 


.27 


20 


.25 


20 


.13 


19 


79 


HETJ1028+3929A 


10 


:28 


43 


.67 


+39: 


29 


:36 


9 


19 


.96 


19 


.90 


19 


.95 


19 


.57 


19 


.49 


HETJ1028+3929B 


10 


:28 


: 14 


.30 


+39: 


:29' 


:34 


.8 


20 


.92 


20 


.81 


20. 


.91 


20 


.60 


20 


.66 


SDSSJ1034+0701A 


10 


:34 


;51 


.47 


+07: 


01 


21 


.2 


19 


.86 


19 


.47 


18. 


.94 


19 


04 


19 


08 


APOJ1034+0701B 


10 


:34 


:51 


.38 


+07: 


:01 


21 





21 


.20 


21 


.26 


21. 


08 


20 


.89 


20 


.99 


2QZJ1056-0059A 


10 


:56 


:44 


.89 


-00: 


:59: 


:33 


,4 


20 


.16 


19 


.92 


19. 


.89 


19 


.80 


19. 


.59 


APOJ1056-0059B 


10 


:56 


:45 


.25 


-00: 


59 


:38 


.1 


20 


.76 


20 


.78 


20 


59 


20 


.58 


20 


33 


SDSSJ1123+0037A 


11 


:23 


:10 


.96 


+00: 


:37: 


45. 


,2 


19 


.02 


19 


.03 


18. 


.90 


18. 


.98 


19. 


.08 


APOJ1123+0037B 


11 


:23 


:07 


,21 


+00: 


:37 


.45. 


7 


20 


.26 


20 


.29 


20. 


.10 


20 


.12 


20. 


.22 


APOJ1225+5644A 


12 


:25 


:45 


.73 


+56: 


:44: 


:40. 


,7 


20 


.05 


19 


.28 


19 


.44 


19 


.35 


19. 


.07 


APOJ1225+5644B 


12 


:25 


:45 


.24 


+56: 


:44: 


:45 


1 


21 


.08 


20 


.52 


20. 


50 


20 


.35 


19. 


.80 


SDSSJ1254+6104A 


12 


54 


;21 


.98 


+61: 


01. 


22. 





19 


.24 


19 


.06 


19. 


01 


18 


92 


18 


.74 


APOJ1254+6104B 


12 


:54 


:20 


.52 


+61: 


:04: 


:36. 





19 


.67 


19 


.56 


19. 


47 


19 


.28 


19. 


.14 


APOJ1259+1241A 


12 


:59 


:55 


.62 


+ 12: 


:41: 


53 


.8 


20. 


.33 


19 


.99 


19. 


93 


19 


.74 


19 


.53 


APOJ1259+1241B 


12 


:59: 


:55 


.46 


+ 12: 


:41 


:51 


.0 


20 


.22 


19 


.90 


19 


87 


19 


.79 


19. 


56 


APOJ1303+5100A 


13 


:03 


:26 


.17 


+51: 


:00' 


:47. 


.5 


20. 


.46 


20 


.33 


20. 


28 


20 


.05 


20. 


.03 


APOJ1303+5100B 


13 


:03 


:26 


.13 


+51: 


00 


:51 


3 


20 


.87 


20 


.60 


20 


66 


20. 


.38 


20 


59 


SDSSJ1310+6208A 


13 


:10 


:37 


89 


+62: 


:08: 


:21. 


6 


18 


.87 


18 


.77 


18. 


63 


18. 


.57 


18 


.35 


APOJ1310+6208B 


13 


:10 


:31 


.96 


+62: 


:08: 


:43. 


.5 


20 


.65 


20 


.49 


20. 


35 


20 


.15 


20 


.12 


SDSSJ1337+6012A 


13 


:37 


:13 


.13 


+60: 


:12: 


06. 


,7 


18 


.68 


18 


.57 


18 


56 


18. 


.34 


18. 


.42 


APOJ1337+6012B 


13 


:37 


:13 


.08 


+60: 


12 


:09 


8 


20 


.15 


20 


.01 


20 


.03 


19 


66 


19 


.80 


SDSSJ1349+1227A 


13 


:49 


:29 


.84 




:27: 


07. 





17 


.92 


17 


.76 


17. 


.71 


17. 


.46 


17. 


.45 


APOJ1349+1227B 


13 


:49: 


:30 


.00 


+ 12: 


:27 


:08 


.8 


19 


.50 


19 


.26 


19 


.10 


18 


74 


18. 


60 


APOJ1400+1232A 


14 


:00 


:52 


.07 


+ 12: 


:32: 


:35 


,2 


20 


.41 


20 


.28 


20. 


.27 


20. 


.13 


19. 


.88 


APOJ1400+1232B 


14 


:00 


:52 


.56 


+ 12: 


:32: 


:48. 





20 


.54 


20 


.42 


20. 


.45 


20. 


.27 


19. 


.99 


SDSSJ1405+4447A 


14 


:05 


:01 


.94 


+44: 


■47 


:59. 


9 


18 


.68 


18 


.14 


17. 


93 


17 


.86 


17. 


.66 


APOJ1405+4447B 


14 


:05 


:02 


.41 


+44: 


:47 


.54. 


,4 


20 


.61 


20 


.12 


19. 


.96 


19 


.89 


19 


60 


APOJ1409+3919A 


14 


:09 


53 


.74 


+39: 


:19 


:60. 


.0 


20 


.31 


20 


.24 


20 


.14 


20 


.06 


19 


.78 


APOJ1409+3919B 


14 


:09 


:53 


.88 


+39: 


:19: 


:53 


,4 


20 


.89 


20 


.92 


20 


87 


20. 


.67 


20 


.31 


APOJ1530+5304A 


15 


:30 


:38 


.56 


+53: 


:04: 


:04. 


2 


20 


.86 


20 


.54 


20. 


42 


20 


.13 


20 


.06 


APOJ1530+5304B 


15 


:30 


:38 


.82 


+53: 


01 


:00 


8 


20 


.67 


20 


.66 


20 


.53 


20 


31 


20 


21 


SDSSJ1546+5134A 


15 


:46 


:10 


.55 


+51: 


■31: 


:29. 


.5 


22 


.34 


20 


.57 


20. 


31 


20. 


.20 


20 


.29 


APOJ1546+5134B 


15 


:46 


:14 


.24 


+51: 


:34 


:05 





20 


.64 


19 


.43 


19 


10 


18 


91 


18. 


.88 


SDSSJ1629+3724A 


16 


:29 


;02 


.59 


+37: 


24 


30 


.8 


19 


.47 


19 


.26 


19. 


10 


19. 


16 


19. 


06 


APOJ1629+3724B 


16 


:29 


:02 


.63 


+37: 


:24: 


:35. 


.2 


19 


.50 


19 


.41 


19 


28 


19. 


.38 


19. 


.31 


SDSSJ1719+2549A 


17 


:19: 


:46 


.66 


+25: 


:49: 


41, 


.2 


20 


.17 


19 


.93 


19 


.90 


19. 


85 


19. 


61 


APOJ1719+2549B 


17 


:19 


:45 


.87 


+25: 


19 


:51 


3 


20 


.08 


19 


.74 


19 


.65 


19. 


65 


19. 


.46 


SDSSJ1723+5904A 


17 


:23 


:17 


.42 


+59: 


04: 


46 


.8 


19 


.40 


19 


.04 


18. 


.95 


18. 


69 


18. 


76 


APOJ1723+5904B 


17 


:23 


:17 


.30 


+59: 


:04: 


43. 


2 


21 


.24 


20 


.43 


20. 


46 


20 


.18 


20. 


.20 


APOJ2128-0617A 


21 


:28 


;57 


.38 


-06: 


:17 


:50 


9 


19 


.84 


19 


.66 


19. 


.87 


19. 


.80 


19 


68 


APOJ2128-0617B 


21 


:28 


:57 


.74 


-06: 


:17 


:57. 


,2 


20 


.03 


19 


.74 


20 


.12 


19. 


.92 


19 


63 


APOJ2214+1326A 


22 


:14 


;27 


.03 


+ 13: 


26 


:57. 





20 


.39 


20 


.19 


20 


.19 


19. 


96 


19. 


60 


APOJ2214+1326B 


22 


:14 


:26 


.79 


+ 13: 


:26: 


:52 


3 


20 


.65 


20. 


.41 


20. 


.26 


19 


.98 


19 


.82 


APOJ2220+1247A 


22 


:20 


:30 


.26 


+ 12: 


:47 


33 


.5 


20 


.11 


20 


.03 


20 


00 


19. 


.88 


19. 


.78 


APOJ2220+1247B 


22 


:20 


:29 


.53 


+ 12: 


:47. 


45 


1 


20 


.92 


20 


.85 


20 


.72 


20 


,34 


20 


21 



A6» 



Ad 



Rr, 



Sample 



3.6 


1.56 


< 200 


22.1 


5.8 


14.1 


2.13 


-1600 


84.5 


17.5 


19.0 


2.30 


-520 


112.4 


10.8 


6.9 


1.64 


< 200 


41.9 


8.8 


18.1 


2.10 


960 


108.4 


5.0 


18.8 


2.12 


-490 


112.4 


14.3 


32.6 


2.53 


290 


189.4 


9.3 


18.6 


0.72 


460 


95.5 


10.8 


3.9 


1.95 


200 


23.8 


18.9 


17.2 


1.51 


< 200 


104.3 


6.5 


22.0 


2.29 


< 200 


129.8 


11.7 


7.5 


1.89 


-1030 


45.6 


6.4 


3.1 


1.25 


850 


18.7 


9.7 


7.2 


2.13 


-350 


43.1 


6.7 


56.3 


1.17 


-330 


332.1 


6.0 


6.0 


2.38 


970 


35.4 


20.7 


17.6 


2.05 


-1010 


105.6 


11.4 


3.6 


2.19 


-840 


21.2 


6.0 


3.8 


1.68 


220 


23.0 


4.0 


46.9 


2.06 


-1850 


281.7 


10.8 


3.1 


1.73 


-610 


18.9 


0.6 


3.0 


1.72 


< 200 


18.3 


12.9 


14.6 


2.05 


1470 


87.7 


0.6 


7.4 


2.23 


1870 


44.2 


3.2 


6.8 


2.08 


480 


40.7 


3.9 


4.1 


1.53 


230 


25.0 


8.8 


42.2 


2.95 


-1450 


236.0 


12.4 


4.4 


0.92 


< 200 


24.5 


15.9 


14.7 


2.17 


-220 


87.5 


6.4 


3.7 


1.60 


-830 


22.7 


3.0 


8.3 


2.07 


-290 


49.7 


15.7 


5.8 


2.00 


-690 


35.2 


19.4 


15.9 


1.99 


1600 


95.5 


43.0 



X 



X' 
photo 

^2 



X' 
photo 

X 2 
X 2 
photo 

„2 



X 

photo 

X 2 
photo 
photo 

■J* 



X 

photo 

X 2 
photo 
photo 

X 2 
X 2 

x 2 
x 2 

photo 
photo 
photo 



Note. — Quasars labeled SDSS or 2QZ are members of the SDSS or 2QZ spectroscopic quasar catalog and are designated as quasar 'A'. 
Quasars discovered from follow up spectroscopy are labeled APO (or HET) and designated 'B'. For pairs discovered from the photometric 
catalog, both quasars are labeled APO (or HET) and 'A' designates the brighter of the two quasars. Extinction corrected SDSS five band 
PSF photometry are given in the columns u, g, r, i, and z. The redshift of quasar 'A' is indicated by column z, Ad is the angular separation 
in arcseconds, Av is the is the velocity of quasar B relative to quasar A in km s , i? pr op is the transverse proper separation in h -1 kpc, 
and x 2 is the value of our color similarity statistic. The last column labeled 'Sample' indicates the selection algorithm used to find the binary. 
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0955+6045 




4000 5000 6000 



7000 8000 9000 10000 
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Fig. 4. — Keck ESI spectra of both members of three binary quasars. The top panel is the binary SDSSJ0955+6045 (z = 0.72, A8 = 18". 6, 
-Rprop = 95.5 h' 1 kpc), the middle panel is SDSSJ1010+0416 (z = 1.51, A6» = 17". 2, i?, pr0 p = 104.3 h' 1 kpc), and the bottom panel is 
SDSSJ1719+2549 (z = 2.17, A6» = 14". 7, i?, pr op = 87.5 h' 1 kpc). The discontinuity in the spectra at 4500A is an artifact of a gap in the 
Echclle orders. 
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Fig. 5. — Keck spectra of two binary quasars. Spectral dissimilarity and the absence of a lensing galaxy in optical and near IR follow-up 
imaging provide strong evidence that these are both binary quasars rather than gravitational lenses. Left: Keck LRIS moderate resolution 
spectra of both members of the binary quasar SDSSJ0048-1051 (z = 1.56, A0 = 3".6, R = 22.1 h" 1 kpc). The absorption feature at 7600A 
is telluric. Right: Keck LRIS low resolution spectra of both members of SDSSJ0248+0009 (z = 1.64, A9 = 6". 8, R = 41.9 h' 1 kpc). 



TABLE 3 

Binary Quasars with A9 < 3" Discovered from Lens Selection 



Name 


RA (2000) 


Dec (2000) 


u 


a 


r 


i 


z 


Ad 


z 


Av 


^?prop 


SDSSJ0740+2926A 


07:40:13.45 


+29:26:48.4 


18.61 


18.46 


18.30 


18.42 


18.48 


2.6 


0.98 


230 


15.0 


APOJ0740+2926B 


07:40:13.43 


+29:26:45.7 


19.98 


19.67 


19.50 


19.68 


20.03 










SDSSJ1035+0752A 


10:35:19.37 


+07:52:58.0 


19.17 


19.13 


18.97 


19.03 


19.14 


2.7 


1.22 


270 


15.8 


APOJ1035+0752B 


10:35:19.23 


+07:52:56.4 


20.62 


20.42 


19.98 


19.84 


19.84 










SDSSJ1124+5710A 


11:24:55.24 


+57:10:57.0 


19.34 


18.66 


18.83 


18.65 


18.44 


2.2 


2.31 


-540 


12.7 


APOJ1124+5710B 


11:24:55.44 


+57:10:58.4 


20.31 


19.83 


19.52 


19.52 


19.42 










SDSSJ1138+6807A 


11:38:09.21 


+68:07:38.8 


18.28 


17.98 


17.87 


17.89 


17.78 


2.6 


0.77 


840 


13.7 


APOJ1138+6807B 


11:38:08.89 


+68:07:36.9 


20.31 


19.74 


19.76 


19.72 


19.58 










SDSSJ1508+3328A 


15:08:42.21 


+33:28:02.6 


17.88 


17.78 


17.81 


17.97 


17.86 


2.9 


0.88 


<200 


16.0 


APOJ1508+3328B 


15:08:42.22 


+33:28:05.5 


20.56 


20.36 


20.15 


20.56 


19.71 










SDSSJ1600+0000At 


16:00:15.50 


+00:00:45.5 


19.23 


19.11 


18.84 


18.95 


19.08 


1.9 


1.01 


-660 


10.6 


NTTJ1600+0000Bt 


16:00:15.59 


+00:00:46.9 






« 21 


»s 21 













Note. — Quasars labeled SDSS are the members of the SDSS spectroscopic quasar catalog and are designated as 
quasar 'A'. Quasars discovered from follow up spectroscopy are labeled APO (or NTT) and designated 'B'. Extinction 
corrected SDSS five band PSF photometry are given in the columns u, o, r, i, and z. These m agnitudes are estimated 
from the deblending algorithm of the main SDSS photometric pipeline IStoughton ~t al. 2002) (except for 1600+0000 
see below). The redshift of the SDSS quasar is given by z, A8 is the angular separation in arcseconds, Av is the the 
velocity of quasar B relative to quasar A in km s _1 ,and R pIO p is the transverse proper separation in h~ 1 kpc. 
t The pair SDSSJ1600+0000A and NTTJ1600+0000B was not deblended by the photometric pipeline because the 
separation is too small. The magnitudes of SDSSJ1600+0000A have contributions from both members of the pair and 
the approximate r and i band magnitudes of NTTJ1600+0000B, were measured from follow up NTT imaging of this 
system. 



plot in Figure [S] The horizontal long-dashed lines indi- 
cate the redshift limits 0.7 < z < 3.0 of the sample, and 
the symbols and dotted and short dashed curves are the 
same as in Figure [3] 

6. CLUSTERING ANALYSIS 

A measurement of quasar clustering on the small scales 
probed by our binary sample 10 hr x kpc < i? pr0 p 
1 h^ 1 Mpc is unprecedented. Our strategy has been to 
overcome the fiber collision limitation and high level of 
shot noise by following up close pairs of quasars to mag- 
nitudes i < 21 fainter than the flux limit of the SDSS 
quasar survey. 

As our pairs have a fainter flux limit than the under- 
lying parent quasar catalog, and the mean density at 
this limit cannot be determined from the brighter sam- 
ple, we cannot use the conventional technique of Monte 
Carlo integration of a random catalog to compute the 
correlation function. However, the mean number den- 



sity of quasar s as a function of magnitude and redshift 
is well known ijCroom et al.ll2004al iRichards et al.ll2005|) 
and can be computed from the quasar luminosity func- 
tion. Below, we explain how we model the quasar lumi- 
nosity function. We then introduce an estimator for the 
quasar correlation function which takes the incomplete- 
ness of our pair survey into account. After determining 
the selection function of our clustering sub-sample, we 
compute its correlation function and compare it to pre- 
vious clustering measurements extrapolated down to the 
scales i?p 10 p ^ 1 h Mpc probed by our binaries. 



6.1. Modeling the Luminosity Function 

At low redshift z < 2.3, the quasar lumin osity func- 
tion has been meas ured by several groups llBovle et alJ 
I2000T iCroom et alJ l2004ri UTichards et all l2005|) . We 
use the double power law B-band luminosity function 




Fig. 6. — SDSS and APO spectra of both members of three binary quasars. The top panel is the binary SDSSJ0245-0113 (z = 2.46, 
AO = 4". 5, .R prop = 26.3 h,- 1 kpc), the middle panel is SDSSJ0740+2926 (z = 0.98, A0 = 2".6, i? pr °p = 15.0 h' 1 kpc), and the bottom 
panel is SDSSJ1124+5710 (z = 2.31, AO = 2". 2, R pTop = 12.7 h' 1 kpc). The binary in the top panel was a member of our overlap sample, 
hence both quasars have SDSS spectra. For both the middle and bottom panels, the black curves are SDSS spectra of the brighter quasar 
in the pair and the red curves are the APO spectra of the fainter companions. The absorption feature at 7600A in the lower two spectra 
is telluric. Although all three of these binaries have separations A9 < 5" characteristic of gravitational lenses, deep optical and near IR 
imaging show no lenses in the foreground. 



(|Bovle et aLl feoOO) 



$(M B , z) 10 o.4(ft+l)[M B -M5(*)] + 10 0.4(/3 fe + l)[M B -M*( 2 )] ' 

(5) 

where j3 x = -1.64, (3 h = -3.43, and $* = 360(/i/0.50) 3 
Gpc -3 mag -1 . The evolution of the break luminosity 



M*(z) follows 

A*5 = M* B (Q)-2.5(k 1 z + k 2 z 2 ), (6) 

with ki = 1.36, k 2 = -0.27, and Af*(0) = -21.15 + 
5 log h. 

The quasar luminosity function is poorly constrained 
between redshifts 2.6 < z < 2.9, because quasar co lors 
cross the stellar locus (see e.g. lRichards et al"ll2001af) . in 
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Fig. 7. — SDSS and APO spectra of both members of three binary quasars. The top panel is the binary SDSSJ1405+44447 (z = 2.23, 
A8 = 7".4, R = 142.8 h" 1 kpc), the middle panel is SDSSJ1530+5304 (z = 1.53, A8 = 4".l, R = 63.3 ft" 1 kpc), and the bottom panel is 
SDSSJ1723+5904 (z = 1.60, A9 = 3". 7, R = 59.0 kpc). The binary in the middle panel is a member of the photometric sample, and 
both spectra were taken at APO. In the other two panels, the black curves are SDSS spectra of the brighter quasar in the pair and the 
red curves are the APO spectra of the fainter companion. Deep imaging of SDSSJ1405+44447 and SDSSJ1723+5904 shows no lens in the 
foreground. The absorption features at 7600A are telluric. 



this range and color selected samples suffer high incom- 
pleteness. As we desire to predict the number density of 
quasars in the redshift range 0.7 < z < 3.0 probed by our 
clustering sub-sample, we devise a simp le interpolation 
scheme to cover the range 2.3 < z < 3.0: iWvithe fc Loehl 
( 20023) use d a simple analytical model to fit the double 
power law lumin osity function in eqn. (J5J to both the 
iFan et all 1)200 1|) high redshift (z > 3.6) luminosity func- 



tion and the iBovle et all l|2000ft low redshift (z < 2.3) 
luminosity functio n. For redshifts z < 2.3 we use the 
IBovle et alJ {2000) expression in eqn. In the range 
2.3 < z < 3.0 we simply linea r ly inter polate between 
eqn. and the IWvithe k. Loebl l(2002a(l fit. Since the 
number of binary quasars in our clustering sample in this 
redshift range is relatively small, our conclusions will be 
insensitive to any uncertainties in this procedure. 
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TABLE 4 

Binary Quasars Discovered in Overlapping Plates and the SDSS+2QZ Catalog"!" 



SDSSJ0012+0052A 
SDSSJ0012+0053B 
SDSSJ0117+0020A 
SDSSJ0117+0021B 
SDSSJ0141+0031A 
SDSSJ0141+0031B 
SDSSJ0245-0113A 
SDSSJ0245-0113B 
SDSSJ0258-0003A 
SDSSJ0258-0003B 
SDSSJ0259+0048A* 
SDSSJ0300+0048B 
SDSSJ0350-0031A 
SDSSJ0350-0032B 
SDSSJ0743+2054A 
SDSSJ0743+2055B 
SDSSJ0747+4318A 
SDSSJ0747+4318B 
SDSSJ0824+2357A 
SDSSJ0824+2357B 
SDSSJ0856+5111A 
SDSSJ0856+5111B 
SDSSJ0909+0002A 
SDSSJ0909+0002B 
SDSSJ0955+0616A 
SDSSJ0955+0617B 
SDSSJ1032+0140A 
2QZJ1032+0139B 
SDSSJ1103+0318A 
SDSSJ1104+0318B 
SDSSJ1107+0033A 
2QZJ1107+0034B 
SDSSJ1116+4118A 
SDSSJ1116+4118B 
SDSSJ1134+0849A 
SDSSJ1134+0849B 
2QZJ1146-0124A 
2QZJ1146-0124B 
SDSSJ1152-0030A 
2QZJ1152-0030B 
SDSSJ1207+0115A 
2QZJ1207+0115B 
2QZJ1217+0006A 
2QZJ1217+0006B 
2QZJ1217+0055A 
2QZJ1217+0055B 
SDSSJ1226-0112A 
2QZJ1226-0113B 
SDSSJ1300-0156A 
2QZJ1300-0157B 
2QZJ1328-0157A 
2QZJ1328-0157B 
2QZJ1354-0108A 
2QZJ1354-0107B 



RA (2000) 


Dec (2000) 


it 


a 


V 


% 


z 


\e 


z 


Aw 


-''•prop 


x 2 


00:12:01.88 


+00:52:59.7 


21.53 


20.82 


20.33 


19.83 


19.63 


16.0 


1.63 


1590 


97.3 


59.8 


00:12:02.35 


+00:53:14.1 


21.03 


20.82 


20.51 


20.22 


20.10 












01:17:58.84 


+00:20:21.5 


17.96 


17.67 


17.82 


17.79 


17.99 


44.5 


0.61 


-260 


212.9 


34.1 


01:17:58.00 


+00:21:04.1 


20.26 


20.01 


20.13 


19.86 


19.89 












01:41:11.63 


+00:31:45.9 


20.23 


20.19 


20.11 


19.87 


19.96 


42.9 


1.89 


1140 


259.2 


2.9 


01:41:10.41 


+00:31:07.1 


20.76 


20.59 


20.50 


20.29 


20.20 












02:45:12.08 


-01:13:14.0 


20.56 


19.85 


19.56 


19.47 


19.33 


4.5 


2.46 


-190 


26.3 


17.0 


02:45:11.90 


-01:13:17.6 


21.14 


20.57 


20.45 


20.38 


20.11 












02:58:15.55 


-00:03:34.2 


18.93 


18.93 


18.66 


18.66 


18.74 


29.4 


1.32 


240 


176.6 


32.9 


02:58:13.66 


-00:03:26.5 


19.67 


19.89 


19.66 


19.75 


19.71 












02:59:59.69 


+00:48:13.7 


19.63 


19.26 


19.22 


19.33 


19.10 


19.6 


0.89 


830 


108.1 


1914.2 


03:00:00.57 


+00:48:28.0 


19.47 


19.01 


16.51 


16.37 


16.05 












03:50:53.29 


-00:31:14.7 


20.35 


20.17 


19.45 


18.99 


18.62 


45.5 


2.00 


-920 


273.8 


492.9 


03:50:53.05 


-00:32:00.1 


19.66 


19.40 


19.32 


19.29 


19.16 












07:43:37.29 


+20:54:37.1 


20.06 


19.86 


19.77 


19.51 


19.43 


35.5 


1.56 


640 


215.5 


15.9 


07:43:36.85 


+20:55:12.1 


20.36 


20.08 


20.04 


19.90 


19.94 












07:47:59.02 


+43:18:05.4 


19.52 


19.24 


19.21 


18.89 


18.75 


9.2 


0.50 


150 


39.9 


6.3 


07:47:59.66 


+43:18:11.5 


19.79 


19.45 


19.36 


19.11 


18.99 












08:24:40.61 


+23:57:09.9 


18.72 


18.51 


18.69 


18.58 


18.59 


14.9 


0.54 


-170 


67.0 


13.2 


08:24:39.83 


+23:57:20.3 


19.00 


18.67 


18.88 


18.70 


18.72 












08:56:25.63 


+51:11:37.4 


18.90 


18.52 


18.64 


18.45 


18.51 


21.8 


0.54 


60 


98.2 


138.9 


08:56:26.71 


+51:11:18.2 


20.03 


19.55 


19.42 


19.17 


19.19 












09:09:24.01 


+00:02:11.0 


16.65 


16.68 


16.61 


16.39 


16.34 


15.0 


1.87 


1700 


90.6 


28.0 


09:09:23.13 


+00:02:04.0 


20.08 


20.06 


20.11 


19.97 


19.82 












09:55:56.38 


+06:16:42.5 


17.79 


18.07 


17.84 


17.81 


17.86 


44.0 


1.28 


-1040 


263.3 


38.6 


09:55:59.03 


+06:17:01.9 


20.35 


20.29 


20.11 


20.21 


20.29 












10:32:44.65 


+01:40:20.5 


18.93 


18.86 


18.84 


18.76 


18.86 


55.1 


1.46 


-1580 


333.5 


16.0 


10:32:43.17 


+01:39:30.0 


20.50 


20.46 


20.27 


20.15 


20.21 












11:03:57.72 


+03:18:08.3 


18.35 


18.36 


18.30 


18.10 


17.96 


57.3 


1.94 


-1730 


345.7 


41.5 


11:04:01.49 


+03:18:17.5 


19.08 


19.04 


19.12 


19.02 


18.97 












11:07:25.70 


+00:33:53.9 


18.98 


18.94 


18.84 


18.51 


18.42 


24.8 


1.88 


280 


150.1 


22.9 


11:07:27.08 


+00:34:07.6 


20.01 


20.04 


20.05 


19.84 


19.80 












11:16:11.74 


+41:18:21.5 


20.35 


18.53 


18.16 


17.94 


17.96 


13.8 


2.99 


890 


76.8 


28.4 


11:16:10.69 


+41:18:14.4 


21.33 


19.44 


19.17 


19.03 


19.03 












11:34:57.74 


+08:49:35.3 


19.30 


19.10 


19.01 


18.83 


18.85 


27.1 


1.53 


-390 


164.4 


6.4 


11:34:59.38 


+08:49:23.3 


19.65 


19.50 


19.34 


19.11 


19.13 












11:46:52.97 


-01:24:46.4 


20.54 


20.23 


20.17 


19.93 


19.65 


28.5 


1.98 


-490 


172.0 


26.0 


11:46:51.19 


-01:24:56.3 


20.53 


20.46 


20.52 


20.35 


20.24 












11:52:40.53 


-00:30:04.3 


18.95 


18.80 


18.93 


18.86 


18.80 


29.3 


0.55 


740 


132.8 


23.1 


11:52:40.10 


-00:30:32.9 


20.32 


20.12 


20.16 


19.99 


20.16 












12:07:00.97 


+01:15:39.4 


19.03 


18.94 


18.84 


18.87 


18.85 


35.4 


0.97 


-260 


200.1 


3.5 


12:07:01.40 


+01:15:04.7 


20.62 


20.40 


20.30 


20.37 


20.36 












12:17:36.18 


+00:06:57.4 


19.96 


19.93 


20.06 


19.76 


19.62 


51.8 


1.78 


-200 


314.0 


13.5 


12:17:35.03 


+00:06:08.6 


19.17 


19.30 


19.32 


19.05 


19.17 












12:17:36.95 


+00:55:22.7 


20.33 


20.24 


20.12 


20.23 


19.84 


40.5 


0.90 


-60 


224.7 


18.3 


12:17:34.25 


+00:55:22.2 


20.25 


19.92 


19.78 


19.82 


19.98 












12:26:24.09 


-01:12:34.5 


17.47 


17.39 


17.23 


17.27 


17.28 


50.6 


0.92 


100 


282.4 


22.8 


12:26:25.58 


-01:13:19.9 


19.74 


19.82 


19.75 


19.79 


19.53 












13:00:45.56 


-01:56:31.8 


18.26 


18.26 


18.14 


17.88 


17.90 


44.5 


1.62 


480 


270.7 


37.6 


13:00:44.52 


-01:57:13.5 


19.81 


19.72 


19.73 


19.59 


19.70 












13:28:30.14 


-01:57:32.8 


19.86 


19.49 


19.52 


19.55 


19.42 


52.6 


2.37 


-890 


309.2 


26.8 


13:28:33.64 


-01:57:27.9 


20.46 


19.92 


19.78 


19.77 


19.75 












13:54:40.40 


-01:08:45.4 


19.48 


19.54 


19.46 


19.29 


19.21 


55.5 


1.99 


740 


334.2 


37.1 


13:54:39.97 


-01:07:50.3 


20.40 


20.13 


19.99 


19.85 


19.64 













Note. — Quasars labeled SDSS or 2QZ are members of the SDSS or 2QZ spectroscopic quasar catalog. The brighter of the two 
quasars is designated 'A', except for SDSS-2QZ pairs for which the SDSS quasar is designated 'A'. Extinction corrected SDSS five 
band PSF photometry are given in the columns u, g, r, i, and z. The redshift of quasar 'A' is indicated by column z, A9 is the 
angular separation in arcseconds, Ad is the velocity of quasar B relative to quasar A in km s _1 , i? pr op is the transverse proper 
separation in h~ 1 kpc, and x 2 is the value of our color similarity statistic. 

t We publish all quasars with proper transverse separations i? pr0 p < 1 h _1 Mpc, however only those with 9 < 60" are listed in 
this table. The entire sample of pairs is published in the electronic version of this article. 

* The quasar SDSSJ0300+0048B has a large BAL trough which explains the very large \ 2 — 1914.2 for this pair. 



Although the luminosity functions we are considering 
are expressed in terms of Mb, both the SDSS spectro- 
scopic survey and our binary quasar survey have flux 
limits in the i-band. Thus in order to compute the 
number of quasars above our flux limit, we need to 



know the cross filter K-correction Kbi(z), between ap- 
parent magnitude i and absolute magnitu de B (see e.g. 
iHogg. Baldrv. Blanton. fc Eisenstem1l2002|) 

M B = i - DM(z) - K Bl (z) (7) 
where DM(z) is the distance modulus. We compute 
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Fig. 8. — Range of redshifts and proper transverse separations 
probed by the binary quasars in our clustering sub-sample. The 
blue circles are binary quasars identified in the SDSS spectroscopic 
sample (8 > 60"). The short dashed line indicates the transverse 
separation corresponding to 9 = 60". The open blue circles indi- 
cate pairs from the overlapping plates. These pairs are required 
to also meet the \ 2 selection criteria (eqn. The magenta cir- 
cles are members of the lens sample and red circles are from \ 2 
sample. The dashed curve indicates the transverse separation cor- 
responding to 8 = 3" below which binaries are found with our 
lens algorithm. The horizontal long dashed curves at z = 0.7 and 
z = 3.0 indicate the redshift limits of the clustering sub-sample. 
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Fig. 9. — Expected mean number of companions per quasar 
with proper transverse separation R prop < 50 h~ l kpc and ve- 
locity difference | A ^j < 2000 km s _1 , as a function of redshift 
calculated from eqn. 1161 which assumes the large scale quasar cor- 
relation function can be extrapolated as a power law to small 
scales. We used the correlation function parameters, 7 = 1.53 
and ro = 4.8 /i" 1 Mpc, measured by PMN from the 2dF quasar 
survey. A perfect survey, with no sources of incompleteness, is as- 
sumed. In a sample of ~ 50, 000 quasars, the predicted number 
of binaries with R < 50 h" 1 kpc is ~ 1. However, our clustering 
sample already contains 20 quasars with transverse separations this 
small, which is compelling evidence for excess clustering over what 
is expected from extrapolating the correlation function. 



TABLE 5 
Summary of Binary Quasar 
Sample 



Algorithm 


Number of Binaries 


Lens 


6 


x 2 


21 


Photometric 


12 


Overlap 


26 


Spectroscopic 


153 



Note. — Number of binary 
quasars with i?, pr op < 1 Mpc, se- 
lected by the various algorithms dis- 
cussed in §|3 and plotted in Figurel3l 



KbAz) from the SPSS c omposite quasar spectrum of 
Vand en Berk et alJ (j2001f) and the Johnson-B and SDSS 
i filter curves. 

The number density of quasars brighter than the flux 
limit i' is an integral over the luminosity function 



n(z, i < i) 



ML 



dM B ^(M B ,z), 



(8) 



where M B brlgM (z) and M B ' (z) are the absolute magni- 
tudes corresponding to the bright and faint end apparent 
magnitude limits as per eqn. J7J. 

To check our model, we compute the cumulative num- 
ber magnitude counts n(< i) for the redshift range 
0.3 < z < 3.0 and compar e to the measurement ove r 
the same redshift interval bv lVanden Berk et al.l {2005). 
We find that our model slightly over estimates the num- 
ber co unts of quasars. Specifically, iVanden Berk et aT] 
(2005) measured n(< 18.5) = 3.74 deg -2 whereas our 
model predicts n(< 18.5) = 4.63 deg -2 . We thus scale 
our model luminosity function down by this ratio so as 
to give the correct number counts for i < 18.5. 

6.2. Estimating the Correlation Function 

If it is true that quasars depart from the gravitational 
clustering hierarchy because of dissipative encounters, 
this would be expected to occur for separations smaller 
than some length scale characteristic of tidal effects or 
mergers. Proper rather than comoving coordinates are 
appropriate for such an investigation. However, clus- 
tering measurements are typically carried out in comov- 
ing coordinates, as these are most intuitive in the linear 
regime where objects are sill moving with the Hubble 
flow. We will compute the correlation function in both 
proper and comoving units, but for definiteness we use 
comoving units in the equations which follow. 

We consider quasars with a maximum velocity differ- 
ence of |Au| < 2000 km s , thus we will measure the 
redshift space correlation function projected over this ve- 
locity interval 



w p (R,z) 



£ S (R, s, z)ds 



(9) 



-77— 7 



where £ s is the quasar correlation function in redshift 
space, w max = 2000 km s , H{z) is the expansion rate 
at redshift z, and the factor of a = 1/(1 + z) converts the 
redshift distance to comoving units. 
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The redshift space correlation function, £ S (R, s,z) is 
the convolution of the velocity distribution in the redshift 
direction, F v (v z ), with the real space correlation function 

£(r,z), 

C(R,8,z)= [ £ (\/R 2 + x 2 ,z) F v (H(z)[x - s})dx. 

(10) 

The radial velocity distribution F v (v z ) has contributions 
from both peculiar velocities and uncertainties in the sys- 
temic redshifts of the quasars. Provided that the distance 
in redshift space over which we project contains most of 
the area under F v , it is a good approximation to replace 
the redshift space correlation function, £ s , under the in- 
tegral in eqn. JJjJ with the real space correlation function, 
£, since radial velocities will simply move pairs of points 
within the volume. 

The small number of close pairs will limit the number 
of bins in R for which we can measure the correlation 
function. Since w p may change significantly over these 
large bins, we choose to measure the volume averaged 
projected correlation function integrated over each ra- 
dial bin. We denote this dimensionless quantity by W p . 
Within a comoving bin [-R m in, -Rmax] it can be written as: 

J^~L Jn™* £(R, a, z)2wRdRds 

VshcW 

where V^heii is the volume of the cylindrical shell in red- 
shift space 

Khell = ^(-Rmax ~ R min) ( ~ jjT^j J ■ ( 12 ) 

The correlation function of a point process is computed 
by comparing the number of pairs detected to the number 
expected in the absence of clustering, taking into account 
the sample limits. We choose the estimator 

1 + W p (R min , R max ) = 49§v ( 13 ) 

Usually, random catalogs are constructed to determine 
the average number of data-random pairs (QR). How- 
ever a subtlety arises in the current context because for 
AO < 60", we have targeted the companion quasars 
(i.e the quasar discovered from follow-up observations) 
to fainter magnitudes than the quasars from the parent 
spectroscopic sample. Rather than estimate the corre- 
lation function from the number of pairs, we use the 
number of companions. Specifically, (QQ) is defined 
to be the number of companions about quasars in the 
parent sample in a given comoving transverse radial bin 
[-Rmin, -Rmax], and (QR) is the average number of random 
companions in this radial bin in the absence of cluster- 
ing. In the case where the parent and companion samples 
are distinct (as is the case for the lens and % 2 samples) 
then (QQ) is just the number of binaries in the bin in 
question. However, if the parent and companion samples 
are identical (as is the case for the overlap and spectro- 
scopic samples), then (QQ) is twice the number of pairs, 
since each of the two quasars in the parent sample have 
a companion. 

Our model of the luminosity function in the previous 
section is used to compute the average number of ran- 
dom companions R about the quasars Q in each parent 



quasar sample, taking into account the flux limits and 
various sources of incompleteness. We separately com- 
pute quasar-random contribution for each selection algo- 
rithm used to define our clustering sample and we then 
add the results. Specifically, for a comoving transverse 
radial bin [i? m in, -Rmax], centered on R, we can write 

JVqso 

(QR) = n{z h i < i')V shell S( Zj , Sj), (14) 

3 

where n(z, i < i') is the number density of quasars above 
the flux limit i' , 14hcii is the volume of each bin in 
eqn. I|12[l. and the sum is over all quasars in the parent 
sample. The quantity S(zj,9j) is the selection probabil- 
ity of detecting a companion about the jth quasar. It is 
a function of angle and redshift since these quantities pa- 
rameterize our various selection algorithms, as discussed 
in § Here 0j = R/D(z) is the angle onto which the 
(logarithmic) center of the bin projects for a quasar at 
redshift z, and D(z) is the comoving distance. 

In what follows we will compare our measured corre- 
lation function to previous larger scale (1 h^ 1 Mpc — 
30 h^ 1 Mpc) clustering measurements extrapolated as a 
power law down to t he scales probed by our binaries. 
Since W p in eqn. I jlip is a function of redshift, we must 
average it over the redshift distribution of our quasar 
sample before a comparison can be made to the mea- 
surement from eqn. (| 1 31) . 

JV qso _ 

Wp(i? m i n , -Rmax) = "TT WpC^nunj -Rmax, Zj) (15) 

PMN found that £(r) is well fit by a power law, with 
7 = 1.53 and r = 4.8 h~ x Mpc, for the 2dF sample taken 
as a whole. They also found ro increases with redshift, 
from ro = 3.4 — 5.9 h~ l Mpc, which we take into account 
in eqn. i|15|) . 

In Figure we show the prediction for the number of 
companions in the bin < i? pl0 p < 50 h~ 1 kpc with 
| At; | < 2000 km s _1 as a function of redshift 

N c = n{z,i < 21)Hheu[l + Wp(Q,50 /T 1 kpc, z)}. (16) 

For only this figure, we ignore the redshift evolution 
of the correlation length and set ro = 4.8 hr 1 Mpc. 
The expected mean number of companions per quasar 
with R < 50 h^ 1 kpc is ~ 2 - 3 x 10~ 5 without in- 
completeness, so that in a sample of ~ 50, 000 quasars 
we expect to find roughly ~ 1 quasar pair with sepa- 
ration < 50 hT 1 kpc. However, our clustering sample 
already contains 20 quasars with transverse separations 
this small. Although our survey is far from complete, 
there is already evidence for excess clustering over what 
is expected from extrapolating the correlation function 
power law. In the next two sections we will make this 
argument more precise. 

6.3. Computing the Selection Function 

To estimate the correlation function, we will compute 
the total number of random pairs expected from eqn. Ijl4|l 
by summing over the redshifts of the quasars in each par- 
ent sample around which we searched around for candi- 
date companions. However, we must first consider sev- 
eral sources of incompleteness. For the x 2 selected sam- 
ple, our x 2 < 20 cut results in a completeness fraction, 
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C(z 7 x 2 < 20), shown in Figure [21 In addition, only a 
fraction of the pair candidates which satisfy the criteria 
for our lens and x 2 selection algorithms have been ob- 
served to date, and the fraction of candidates observed 
varies with angular separation, because of our tendency 
to observe small separation pairs first. Finally, for the 
9 > 60" pairs which we find in the SDSS spectroscopic 
sample, the completeness fraction is just that of the SDSS 
quasar survey and does not vary with angle. 

The selection probability in eqn. (|14|) can therefore be 
written 



Fleas (#) 



< 3" 



S(z, 8) = \ F X 2 {6)C (z\ x 2 < 20) 3" < 6 < 60" (17) 



6 > 60" 



where F\ eus {6), F x 2(0), and F spcc are the completeness 
fractions of the lens, x 2 > an d spectroscopic samples, re- 
spectively. 

We take F spcc = 0.95 following IVanden Berk et all 
(2005). The angular selection functions, Fi ens (9) and 
F x 2(9), of the lens and x 2 selection algorithms, are com- 
puted by comparing, in bins, the number of pair can- 
didates which have been observed to date to the total 
number of candidates. 

For 8 < 15", we choose the bin spacing such that each 
bin contains at least 6 objects. For 15" < 9 < 60" we 
use ten logarithmically spaced angular bins. 

The selection probability for the angular bin [9k, &k+i] 



is 



F, 



N 



bscrvcd 



(18) 



remaining 



There is an uncertainty in the number of remaining 
candidates for the lens selection, which is based on the 
goodness of fit of a m ulti-component PSF to the im- 
ages of SDSS quasars l|Pindor et al.ll2003t llnada et al.1 
2003). Because our observations were heavily biased to- 
wards those candidates which were likely to be quasar 
pairs, we bracket the angular selection function F\ cns (6) 
with upper and lower limits corresponding to pessimistic 
and optimistic assessments of the number of remaining 
candidates which are likely to be confirmed as binaries. 
The lens candidates are given a grade of A, B, or C de- 
termined by the goodness of fit of the multi-component 
PSF. For the lower limit on Fi ens (9), N Tema i n i ng is taken to 
be all candidates, whereas for the upper limit -/V rema i n i ng 
is restricted to only the candidates which received a grade 
of A. Note that to be conservative, we have ignored the 
intrinsic incompleteness of the lens selection algorithm, 
which is a function of se paration, magnitude, and flux 
ratio ijPindor et al J 1200^1 of the pair, and focus only on 
the incompleteness of our observations of the candidates. 

The color-color diagrams of the candidates identified 
by x 2 selection, were visually inspected prior to observa- 
tions candidates which overlapped the stellar locus were 
given a lower priority. However, some interlopers were 
observed depending on a variety of criteria. For instance, 
we were more likely to observe pairs with particularly 
small angular separations, or redshifts z > 2 (because of 
their use for Lya forest studies), or those with particu- 
larly bright magnitudes. Including the stellar locus inter- 
lopers identified by the x 2 statistic as part of iV rema i n i ng 
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Fig. 10. — Angular selection function of binary quasars with 
8 < 60". We have pasted F lens (0) and F x2 (0) together at 8 = 3". 
For 8 < 3", the red and blue histograms represent, the upper and 
lower limits on the fraction of lens pair candidates observed to 
date, respectively. The green histogram shows the fraction of \ 2 
pair candidates observed so far. 



in eqn. (|18|l would overestimate the number of candidates 
likely to be quasar pairs and hence underestimate our se- 
lection function. To this end, we apply additional criteria 
to these remaining candidates to filter out color-matches 
which arc likely to be quasar-star pairs. However, all 
quasar-star pairs we confirmed spectroscopically are in- 
cluded as part of iV bserved in eqn. Ijl8|l . In this way, we 
conservatively err on the side of overestimating our se- 
lection completeness, or underestimating the correlation 
function. 

In addition to the criteria in eqn. ifKjl. we only in- 
clude x 2 candidates which either (A) have colors con- 
sistent with the UV-excess region of color space (it — g < 
0.6) but outside the r egion populated by white dwarfs 
l)Richards et al.ir2002a|L (B) matched a member o f the 
faint photometric catalog of iRichards et all l)2004|) (ex- 
tended to the DR3 region), or (C) were optically unre- 
solved and matched a FIRST radio source. Objects that 
satisfy these criteria have a > 90% probability of being a 
quasar. For the UV-excess criteria, we conservatively re- 
quired that the candidates' la photometric errors left it 
inside the UV-excess region and outside the white dwarf 
region. 

The angular selection function of our binary quasar 
survey is shown in Figure ITU1 

6.4. Excess Clustering 

The individual quasar-random contributions for each 
selection algorithm were computed by carrying out the 
sum in eqn. i|14|) over the respective parent samples. 
Specifically, we summed over 39, 142 quasars which were 
the parent sample for the lens sample, the 59, 608 
(SDSS+2QZ quasars) which were the parent sample of 
the x 2 algorithm, and the 52,279 (SDSS only) quasars 
which made up the parent sample of the spectroscopic 
sample. 

For the proper (comoving) projected correlation func- 
tions, we binned our clustering sample into 15 logarith- 
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R (kpc/h) R (kpc/h) 

Fig. 11 , — Small scale quasar clustering in proper coordinates Left: Comparison of the projected correlation function W p (R m i n , Rmax, z) 
(see eqn. llll measured from our clustering sample with the prediction of the large scale measurement of PMN, extrapolated as a power law 
down to the scales probed by our binaries. Error bars are one sigma Poisson counting errors. The blue points show the measurement of 
the pr ojec ted correlation function if we use the lower limit (optimistic) for the selection function of our lens algorithm (blue histogram in 
Figure lTOl to predict the number of quasar-random pairs, rather than the upper limit, which are shown by the black points. The blue points 
are offset slightly to the right for the sake of illustration. Red rectangles indicate the prediction based on the large scale measurement of 
PMN, where the height of the rectangles indicate the range of predictions based on one sigma errors in the correlation length measurements, 
and the width indicate the bin used for each measurement. Right Ratio of the projected correlation function to the (best-fit) prediction of 
PMN. 
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Fig. 12. — Same as Figure UTI except in comoving units. 



mically spaced bins in the range 9 h 1 kpc < -R pr0 p < 
3 hr 1 Mpc (20 hr 1 kpc < R < 7 hr 1 Mpc). Both pro- 
jected correlation functions are compared to an extrap- 
olation of the large scale power law measured by PMN 
in the left panels of Figures El and El The PMN cor- 
relation function is in comoving coordinates, which was 
translated to proper coordinates in each element of the 
average in eqn. (|15f) . The error bars are one sigma Pois- 
so n counting erro rs, where we used the fitting formula 
in iGehrelsl l)1986|) for N < 30. The blue points show 
the measurement of the projected correlation function 
if we use the lower limit (optimistic) for the selection 
function of our lens algorithm (blue histogram in Fig- 
ure I10|) to predict the number of quasar-random pairs, 
while the black points use the upper limit. Uncertainty 
in our selection function changes only the innermost few 
bins, because our lens algorithm is restricted to angular 
separations 9 < 3". These blue points are offset slightly 



to the right for the sake of illustration. The red rectan- 
gles indicate the prediction for the projected correlation 
function based on the large scale measurement of PMN, 
extrapolated as a power law down to the scales probed by 
our sample. The heights of the rectangles are the range 
of predictions based on the one sigma errors on the cor- 
relation lengths (for each redshift interval) published by 
PMN, and the widths of the rectangles indicate the bin 
used for each measurement. The ratio of the measured 
projected correlation functions to the extrapolation of 
the best fit prediction of PMN are shown in the right 
panels of Figures ^2 an d El 

A summary of the data plotted in these Figures is given 
in Tables and In the Tables we compare our calcu- 
lation for the number of random pairs in each bin (QR) 
(eqn. I14|) for our survey to (Qi?) pe rfoct, which is the ex- 
pected number of random pairs for a 'perfect' survey with 
no sources of incompleteness. A comparison of these two 
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quantities is a measure of our incompleteness in each bin. 

Our clustering measurement agrees with the measure- 
ments of PMN for large proper (comoving) separations 
Rprop > 400 h- 1 kpc (R > 1 h" 1 Mpc), where the PMN 
measurement is valid. However, we detect a significant 
excess over the expectation from extrapolating the power 
law of the correlation function to smaller scales. Al- 
though the uncertainties are large because of shot noise 
(and our uncertain selection function for 9 < 3"), the 
excess is an order of magnitude for proper (comoving) 
separations i? pr0 p < 40 h^ 1 kpc (R < 100 h^ 1 kpc). 
The excess is largest (~ 30) on the smallest scales 
-Rprop ~ 15 h- 1 kpc (R ~ 30 br 1 kpc). 

Our clustering statistic W p (R m i n , R max ) is the average 

over a cylinder shell with velocity extent 4000 km s~ , 
corresponding to a typical extent r prop ~ 18 h^ 1 Mpc 
(r ~ r )~ 7 , W p (R min , R max ) progressively flattens to- 
ward small scales, as is seen in red rectangles in the left 
panels of Figures ll II and IT2*1 In contrast, our measured 
Wp(R m i n , Rmax) does not flatten, but rather resembles a 
power law in projection. This necessarily implies that the 
true three dimensional correlation function of quasars is 
not a power law £(r), but rather gets progressively steeper 
on small scales. 

7. DISCUSSION 

On large proper comoving scales > 1 h^ 1 Mpc, quasars 
at z ~ 1.5, similar to high redshift galaxies, are strongly 
biased, and have nearly the same correlation length that 
galaxies do in the local universe. In the previous sec- 
tion we argued that the correlation function of quasars 
becomes significantly steeper on sub-Mpc scales. Does 
small scale galaxy clustering or (quasar-galaxy cluster- 
ing) show a similar trend? As we compare to previ- 
ous clustering measurements in this section, all distances 
quoted are comov ing. 

At low redshift. iZehavi et alJ <|2004f) measured the pro- 
jected correlation function of the nearby (z < 0.2) SDSS 
main galaxy sample (L ~ £*), down t o scal es as small 
as ~ 100 h^ 1 kpc. lEisenstein et al.l l)2005|) measured 
the small scale clustering of Luminous Red Galaxies 
(L ~ 5 - 10L*) at intermediate redshift (0.2 < z < 0.3), 
down to scales of ~ 300 hr 1 kpc, by cross correlating 
them with L* galaxies. Both of these studies find small 
deviations from a power law 7 ~ —1.9, and the correla- 
tion function does steepen toward smaller scales. How- 
ever, if these correlation function power laws were ex- 
trapolated from large (~ 1 h^ 1 Mpc) scales to small, 
as we have done above for quasars, the clustering 'ex- 
cess' would not be more than ~ 30%. In comparison, we 
have found a much larger clustering excess, of a factor 
of ~ 7 on scales R ~ 100 ft. -1 kpc and ~ 2 — 3 on scales 
R ~ 300 h^ 1 kpc (see Table[7||. One could argue that the 
proper compariso n is with the square of th e clustering ex- 
cess measured bv lEisenstein et all l)2005|) (see discussion 
below) , since it is a cross-correlation with less luminous 
galaxies. However IZehavi et alJ i|200ft) measured the au- 
tocorrelation of LRGs at ~ 400 ft -1 kpc, and also found 
an excess smaller than ~ 30%. 

Another complementary set of low redshift observa- 
tions which probe clustering on small scales (lOft -1 kpc— 
l/t" 1 Mpc) are the SPSS galaxy-galaxy le nsing studies 
l|Guzik fc Seliakll2002l iSheldon et al.ll200l|) . which mea- 
sure the galaxy-dark matter cross-correlation function. 



No excess small scale clustering over a power law was 
detected in these studies either. 

A comparison with small scale galaxy clustering at 
high redshift is clearly mor e relevant to our purposes. 
lAdelberger et alJ l|2003fl and lGoil et all (|2004ft have mea- 
sured the correlation function of galaxies at z ~ 3 
and z ~ 1, respectively, down to scales as small as 
~ 100 — 200ft. -1 kpc, and no small scale excess cluster- 
ing was found. As the redshift ranges of these measure- 
ments bracket the redshift range probed by our sample 
of binary quasars, we expect that galaxies in the redshift 
desert z ~ 1.4 — 2.5, coeval with the bulk of our binary 
quasars (see Figure |2J , also do not show any enhance- 
ment in small scale clustering^ 

Thus, as first suggested bv lDiorgovskil l)1991[K quasars 
at z ~ 1.5 depart from the power law clustering hi- 
erarchy followed by galaxies, both in the local uni- 
verse and at high redshift. Galaxy interactions are 
often implicated as a means of triggering, fueling, or 
forming active galactic nuclei iTT bomre fc Toomrelll972t 
iBarnes fe Hernouistl Il996t iBahcall et al.l Il997j) . and it 
has been claimed by several authors that excess quasar 
clusterin g is the hallmark o f these dissipative interaction 
event s llDiorgov skil Il9 91t iKoc hanek. Falco. fc Munod 
ll999HMortlock. Webster, fc Francis! 119991) . However, if 
this is the similar small scale enhancement is to 

be expected in the quasar-galaxy correlation function, as 
follows. 

If the quasar-quasar correlation function is £qq, and 
quasars trace galaxy overdensities with a linear relative 
bias 6q defined by Sq = bc^Sa, then the ratio of the 
quasar-quasar correlation to the galaxy-galaxy correla- 
tion is £qq/£gg = Oq, and the ratio of the quasar-galaxy 
cross correlation function to galaxy-galaxy clustering is 
Cqg/Cgg = bQ. Since we have argued that small scale 
quasar clustering is enhanced relative to galaxy cluster- 
ing at the same redshift, we expect to see the square root 
of that enhancement in the quasar-galaxy cross correla- 
tion. 

It has long been known that quasars are associated 
with enhancements in the dist r ibution of galaxies 
(iBahcall. Schmidt, fc GiT^ml Il969t lYee fc Greer] JlQSl 
| 1987tlBahcall fc Chokshill99lHLaurikainen fc Saloll99ft 
I Smith. Bovle. fc Maddoxl 1199ft IHall. Green, fc Coheri 
[Tgga IHa.ll fc Oreenl 1199ft ISmitn. Bovle. fc Maddo* 
2000), and the foregoing argument has led us to ask 
if there is a small scale excess of galaxy pairs around 
quasars. However, nearly all studies of quasar-galaxy 
correlations have been restricted to low redshifts 
(z < 0.6) and therefore low-luminosity, and the only cor- 
relation function measurements that probed the redshift 
range 1 < z < 3 relevant for comparison to our sample of 
binaries measur ed a marginal cross-correlation of qu asars 
with galaxies llTeplitz. McLean, fc MalkarJ I1999D. or 
none at all ilBovle fc Couchl 119931: lOroom fc Shanks! 
1999; Infa nte et alJ 120031) . We note that the disagree- 
ment in the literature on whether radio-loud quasars 
are located in richer environments llYee fc Green ] |1984|. 
19871 lEllingson. Yee. fc Greenlll99l or not llFisher et al. 
1996; jMcLure fc Dunlonl 120011 iFinn. Impev. fc Hooper 
2001 ; iWold et al.ll2001[) is not an issue in the current 
context, since only only one close binary in our clustering 
sub-sample, LBQS1429— 008, has a radio loud member. 
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TABLE 6 

Projected Correlation Function in Proper Coordinates 



R ■ 


flmax 


{QQ) 


{QR) 


(QR) erfect 
per cc 


W p 


w PMN 

P 


Ratio 


9.00 


13.26 


4 


0.02014 (0.01136) 


0.03147 


197.60 (350.98) 


7.89 


25.0 (44. 


13.26 


19.53 


7 


0.06694 (0.05749) 


0.06922 


103.57 (120.76) 


6.40 


16.2 (18. 


19.53 


28.76 


5 


0.04506 


0.2240 


109.96 


5.17 


21.3 


28.76 


42.37 


2 


0.09635 


0.4860 


19.76 


4.18 


4.73 


42.37 


62.40 


3 


0.11784 


1.055 


24.46 


3.37 


7.26 


62.40 


91.92 


2 


0.1654 


2.288 


11.09 


2.71 


4.09 


91.92 


135.39 


6 


0.6090 


4.964 


8.85 


2.17 


4.07 


135.39 


199.42 


7 


0.5938 


10.77 


10.79 


1.74 


6.22 


199.42 


293.74 


8 


1.142 


23.37 


6.01 


1.38 


4.36 


293.74 


432.67 


20 


3.305 


37.30 


5.05 


1.09 


4.64 


432.67 


637.31 


23 


12.59 


13.25 


0.83 


0.85 


0.97 


637.31 


938.74 


49 


27.31 


28.75 


0.79 


0.66 


1.20 


938.74 


1382.73 


71 


59.26 


62.38 


0.20 


0.50 


0.39 


1382.73 


2036.71 


195 


128.6 


135.3 


0.52 


0.38 


1.37 


2036.71 


3000.00 


316 


279.0 


293.6 


0.13 


0.27 


0.48 



.5) 

-9) 



Note. — Data for clustering measurements shown in Figure ITTI The 15 logarithmically spaced 
bins of proper transverse separation are given by -R m i n and iJ max in kpc. The number of observed 
pairs in our clustering sub-sample for each bin is given by {QQ). Our calculation of the number of 
random pairs in each bin is {QR) (ean. ITU , and (Q-R) per f G ct is the expected number of random pairs 
for a 'perfect' survey with no sources of incompleteness._ Our measurement of the projected correlation 
function is given by Wp(R m i n , Rmax). The quantity Wp MN (-R m m, -Rmax) is the prediction from the 
larger scale clustering measurements of PMN, extrapolated as a power law down to the scale pr obed 
by our binaries, and averaged over the redshift distribution of our parent quasar sample (see can. 1151 . 
The ratio of our measurement to the prediction from PMN is given in the last column. The quantities 
in parentheses are the measurements if w e use the lower limit for the selection function of the 'lens' 
algorit hm (blue histogram in Figure ITU1 to predict {QR), which corresponds to the blue points in 
Figure ITTI 



TABLE 7 

Projected Correlation Function in Comoving Coordinates 



-^-min 


-Rmax 


{QQ) 


{QR) 


{QR)perfect 


Wp 


" p 


Ratio 


20.00 


29.56 


4 


0.01251 (0.007935) 


0.02452 


318.69 (503.12) 


8.73 


36.5 (57.6) 


29.56 


43.68 


5 


0.03213 (0.02595) 


0.06007 


154.61 (191.66) 


7.06 


21.9 (27.1) 


43.68 


64.54 


5 


0.05895 (0.05623) 


0.1674 


83.81 (87.92) 


5.71 


14.7 (15.4) 


64.54 


95.38 


3 


0.07226 


0.3886 


40.52 


4.61 


8.80 


95.38 


140.95 


3 


0.1127 


0.8487 


25.63 


3.71 


6.91 


140.95 


208.28 


2 


0.1907 


1.853 


9.49 


2.98 


3.18 


208.28 


307.80 


3 


0.4731 


4.047 


5.34 


2.39 


2.24 


307.80 


454.85 


10 


0.6467 


8.839 


14.46 


1.91 


7.58 


454.85 


672.16 


5 


1.032 


18.75 


3.85 


1.52 


2.54 


672.16 


993.29 


21 


3.160 


29.00 


5.65 


1.20 


4.71 


993.29 


1467.84 


20 


10.80 


13.29 


0.85 


0.94 


0.91 


1467.84 


2169.12 


12 


23.71 


24.96 


0.77 


0.73 


1.06 


2169.12 


3205.45 


70 


51.78 


54.51 


0.35 


0.56 


0.63 


3205.45 


4736.89 


166 


113.1 


119.0 


0.47 


0.42 


1.11 


4736.89 


7000.00 


317 


246.9 


259.9 


0.28 


0.31 


0.92 



Note. — Same as Table 151 except in comoving coordinates (Figurc ll2l , 



Historically, there has been significant scatter in the 
low redshift measurements of quasa r-galaxy correlations 
(see|B rown, Bovle. fc Websterl200H Table 1 for a compi- 
lation of recent studies) caused by heterogeneous quasar 
samples, methodology, and imaging depths. Further- 
more, the measurement of a factor of ~ 4 enhancement of 
the number of companion galaxies around quasars com- 
pared to the mean number around galaxies, measured by 
the Hubble Spac e Teles cope Studies of ~ 40 nearby AGN 
by iFisher et all l|1996ft and IMcLure fc Dunlod (l2001lj 
has be en called into question bv lFinn"Tmpe^feHoTperl 
<|2001l> as being the result of sample biases. We thus fo- 
cus on the most recent determinations of quasar galaxy 



correlations, albeit at low redshifts z < 0.3, by the 2dF 
and SDSS; these surveys have samples of ~ 10, 000 AGN 
s urrounded by ~ 100, 00 g alaxies at their dis posal. 

iCroom et~al] l)2004a|) and I Wake et all l|2004|) measured 
the ratio £qg/£gg = &q, in the 2dF and SDSS surveys 
respectively, finding it to b e consistent with un ity on 
small and large scales. The ICroom etaTJ (|2004a) mea- 
surement extends down to ~ 900 h~ x kpc; whereas the 
I Wake et alJ 1)200 4*1 measurement probes down to scales ~ 
400ft- -1 kpc. Both of these studies cross-correlated spec- 
troscopic samples of galaxies with spectroscopic samples 
of AGN; hence, they are limited to low redshift z < 0.3, 
because the spectroscopic galaxy samples are shallow, 
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and large scales, because of fiber collisions. At z ~ 0.2 
the fiber collision limit corresponds to ~ f60 h^ 1 kpc 
in the SDSS and ~ 80 hr 1 kpc in the 2dF. What is 
needed is a clustering study using photometric galaxies 
around spectroscopic quasars which could extend to red- 
shift z ~ 0.5 and resolve the scales R ~ 10 h^ 1 kpc of 
most interest to u s for the excess quasar clustering (see 
iSerber et al.l l2004. in preparation). 

8. SUMMARY AND CONCLUSIONS 

We have presented a sample of 218 new quasar pairs 
with transverse separations i? pr op < 1 h Mpc over the 
redshift range 0.5 < z < 3.0. Of these, 65 have angu- 
lar separations 9 < 60" below the SDSS fiber collision 
scale. Our 26 new pairs with proper separations -R pr0 p < 
50 h^ 1 kpc (9 < 10") more than doubles the number 
of binaries known with splittings this small. Although 
these binaries were discovered with a variety of selection 
algorithms, we defined a statistical sub-sample selected 
with homogeneous criteria, and computed its selection 
function taking into account sources of incompleteness. 
We presented the first measurement of the quasar cor- 
relation function on the small proper (comoving) scales 
10 h- 1 kpc - 400 hr 1 kpc (20 hr 1 kpc - 1 hr^ Mpc). 
We detect an order of magnitude excess clustering for 
proper (comoving) separations i? pr0 p < 40 ft. -1 kpc 
(R < 100 hr 1 kpc), which grows to ~ 30 on the small- 
est scale probed by our sample, i? pr0 p ~ 15 h^ 1 kpc 
(R ~ 30 hr 1 kpc). 

We reviewed recent small scale measurements of galaxy 
clustering and quasar-galaxy clustering and discussed the 
results in relation to the excess small scale quasar clus- 
tering that we measured. The quasar-galaxy correlation 
function of redshift z ~ 1.5 quasars should show a small 
scale clustering enhancement with amplitude roughly the 
square root of the enhancement detected here. How- 
ever, existing studies of the environments of quasars at 
z oj 1.5 have been plagued by small sample sizes and 
lack the statistics to address the clustering strength on 
the ~ 100 h^ 1 kpc scales of interest. 

Deep imaging of the binaries published here will pro- 
vide valuable information about their environments. The 
detection of significant overdensities of galaxies coeval 
with the binary would support the idea that enhanced 
quasar activity is triggered by galaxy interactions, and it 
might suggest that quasars at high redshift trace the bi- 
ased peaks whi ch are the progenitors of the rich clusters 
we see today iRfatathioii fc Rees)l1988t ICole fc KaiseH 
19891 iNusser fc Silkl 1199.1 IDiorgovski et alJ 1199ft 
Diorgovski 1199ft IDiorgovski et all 12003(1 . Intriguinglv. 
Fuku gita et alJ l|2004fl took dee p Subaru image s of th e 
z = 4.25 quasar discovered bv iSchneider et all (|200Di) . 
and found no evidence for an overdensity of galaxies. 
Similar deep imaging studies of high redshift quasar 
en vironments (fo r singl e quasars) have been conducted 
by llnfante et al.l (|2P03j) at z ~ 3, and by Djorgovski 
and collaborators at z > 4 j Diorgovski et al J 1199ft 
I Diorgovski 119991: IDiorgovski et alJ 120031 iStiavelli etall 
2005). The binaries in our sample offer an opportunity 
to conduct analogous studies over a range of lower and 
more accessible redshifts, with the added bonus that one 
expects these extremely rare binary systems to trace 
even richer environments. 

Measurements of the shape of the quasar-quasar 



and quasar-galaxy correlation function on the small 
scales probed by our binaries will yield valuable in- 
sights into the physical processes that trigger quasar 
activity and will help explain how quasars are em- 
bedded in the structure formation hierarchy. Repro- 
ducin g the excess clustering with semi-analytical mod- 
els l|Kauffmann fe Haehneltil2002|) and halo models (see 
e.g. PMN) would provide constraints on the distribu- 
tion of quasars in dark matter halos. Another inter- 
esting question is whether the quasars in these bina- 
ries have significantly longer lifetimes than 'field' quasars 
l|Haiman fc HiJl20?m iMartini fc Weinbergll2001(l . which 
would have interesting i mplications for the ma sses of 
supermassive black holes llWvithe fc Loebll2004|) in the 
richest regions of the Universe. 

We close with the reminder that our survey for quasar 
pairs is ongoing and less than 50% complete. We expect 
to find a comparable number of binaries in the current 
SDSS quasar sample. Furthermor e, the faint pho t omet- 
ric quasar selection techniques of [Richards et all ( 2004) 
aim to construct a sample of ~ 10 6 quasars. Extrapolat- 
ing from the number of pairs published here, we would 
expect to find ~ 1000 new binaries in a sample of this 
size, which would allow a much more precise measure- 
ment of the correlation function on small scales. 
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TABLE 8 

Projected Pairs of Quasars Discovered in Overlapping Plates 



Name 


RA 


Dec 


a 


9 


r 


i 


2 


z 


A6» 


^prop 


SDSSJ0000+0055A 


00:00:42.91 


+00:55:39.5 


18.30 


18.16 


17.99 


18.01 


17.88 


0.95 


170.7 


961.9 


SDSSJ0000+0055B 


00:00:42.91 


+00:55:39.5 


20.01 


19.69 


19.44 


19.39 


19.43 


1.18 






SDSSJ0004+0000A 


00:04:42.18 


+00:00:23.4 


19.24 


19.09 


18.85 


18.95 


19.07 


1.01 


134.0 


626.4 


SDSSJ0004-0001B 


00:04:42.18 


+00:00:23.4 


20.91 


20.39 


20.22 


19.83 


19.52 


0.58 






Note. — Quasars 


labeled SDSS 


or 2QZ are 


members 


of the 


SDSS 


or 2QZ spectroscopic quasar catalog. 



The redshift of each quasar is indicated by the column z and the foreground quasar is always designated 'A'. 
The column labeled R pTO p is the transverse proper separation at the foreground quasar in h —1 kpc, extinction 
corrected SDSS five band PSF photometry is given in the columns u, g, r, i, and z, and A8 is the angular 
separation in arcseconds 



TABLE 9 

Projected Pairs of Quasars Discovered From Follow Up Spectroscopy 



Name 




RA 


Dec 


u 


9 


r 


i 


z 


z 


A9 


R 


x 2 


APOJ0002- 


-0053A 


00:02:12.53 


-00:53:11.7 


20.72 


20.59 


20.34 


20.12 


20.17 


1.54 


7.3 


44.5 


20.9 


SDSSJ0002- 


-0053B 


00:02:12.53 


-00:53:11.7 


20.10 


19.64 


19.52 


19.41 


19.23 


2.21 








SDSSJ0036- 


1109A 


00:36:49.63 


-11:09:29.8 


19.08 


18.85 


18.67 


18.45 


18.67 


1.51 


4.7 


28.5 


19.0 


APOJ0036- 


-1109B 


00:36:49.63 


-11:09:29.8 


20.92 


20.64 


20.46 


20.52 


20.24 


2.18 









Note. — Quasars labeled SDSS or 2QZ are members of the SDSS or 2QZ spectroscopic quasar catalog. Quasars 
discovered from follow up spectroscopy are labeled APO. The redshift of each quasar is indicated by the column z and 
the foreground quasar is always designated 'A'. The column labeled i? pr0 p is the transverse proper separation at the 
foreground quasar in h _1 kpc, extinction corrected SDSS five band PSF photometry is given in the columns u, g, r, i, 
and z, A9 is the angular separation in arcseconds, and \ 2 is the value of our color similarity statistic. 
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APPENDIX 

In this appendix we present the results of all of our follow up observations not included above, as well as list of 
projected pairs of quasars in the SDSS quasar sample which have similar redshifts. These tables will facilitate future 
studies of close quasar pairs and prevent the duplicate observations of candidates already observed by the SDSS or 
our follow up observation program. 

In Table 00 we list all projected pairs of quasars in the combined SDSS+2QZ spectroscopic catalog with proper 
transverse separations i? pr0 p < 1 h^ 1 Mpc at the foreground quasar and redshift difference Az < 0.5. All projected 
pairs of quasars discovered from our follow up spectroscopic observations are listed in Table EJ Finally, Table ITUl lists 
all of the quasar-star pairs identified from our follow up observations. The full versions of these tables are published 
in in the electronic version of this article. 
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RA 
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u 


a 


r 


i 


z 


z 


A8 
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